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PREFACE 


i 


This  report  Is  submitted  to  the  U.S.  Coast  Guard,  Department  of 
Transportation,  by  Flow  Research  Company  (FRC)  in  completion  of  the 
research  work  on  "A  Study  to  Conduct  Experiments  Concerning  Turbulent 
Dispersion  of  Oil  Slicks,"  under  Contract  No.  DOT-CG-G1688-A.  The 
objective  of  the  research  study  was  to  characterize  (qualitatively  and 
quantitatively)  and  to  identify  the  mechanisms  behind  turbulence  due  to 
wind,  waves  and  currents  operating  independently  and  coupled,  with  and 
without  an  oil  slick  present. 

Thi9  research  study  is  a part  of  the  integrated  research  program 
sponsored  by  tjie  Deepwater  Ports  program  and  the  U.S.  Coast  Guard  to 
identify  the  sea-state  threshold  beyond  which  containment  and  recovery 
of  oil  are  impractical  because  of  the  natural  dispersion  of  oil  slicks. 
Other  parts  of  the  research  program  include  (i)  a preliminary  theoretical 
study  by  Arthur  D.  Little  Jnc.  under  Contract  No.  DOT-CG-61-505A  (Raj, 
1977),  (ill  a laboratory  study  by  Massachusetts  Institute  of  Technology 
(MIT)  under  Contract  No.  DOT-CG-61-802A  to  investigate  the  effects  of 
oil  properties  on  oil  dispersion,  and  (iii)  a field  study  by  the  Naval 
Underwater  Systems  Center  to  conduct  observations  of  surface  layer 
turbulence  associated  with  wind  stress  and  breaking  waves. 

The  research  team  at  FRC  was  headed  by  Dr,  J.  T.  Lin  and  assisted 
by  Drs.  M.  Gad-el-Hak  and  H.  T.  Liu.  Mr.  R.  Srnsky's  assistance  in 
conducting  experiments  is  gratefully  acknowledged.  Technical  Services, 
provided  by  M.  Wallace,  D.  Turner,  I.  Smith,  M.  A.  Kerstein  and  M. 

McKenna  in  drafting,  typing  and  editing  is  recognized  gratefully.  We 
also  acknowledge  Drs.  S.  Veenhuizen  and  M.  Rizk  for  their  valuable 
comments  on  the  write-up  of  this  report.  The  research  study  was  moni- 
tored by  Mr.  R.  A.  Griffiths  of  the  Coast  Guard.  We  gratefully  acknowledge 
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1.  INTRODUCTION 

Because  of  the  increasing  demand  for  oil  imports  to  the  United 
States,  transport  facilities  capable  of  unloading  from  the  200,000  to 
700,000  DWT  class  of  very  large  crude  carriers  (VLCC)  ore  urgently 
required  (Patton  and  Johnson,  1974).  A direct  shipment  by  a VLCC  to  a 
coastline,  for  example,  the  Gulf  Coast,  has  considerable  economic 
advantage.  Unfortunately,  the  Gulf  Coast  has  a gradually  sloping  con- 
tinental shelf,  which  requires  that  deep-water  ports  for  unloading  oil 
directly  from  the  VLCC  must  be  located  offshore.  Severe  sea  and  weather 
conditions  in  this  area  have  been  reported.  According  to  the  National 
Climatic  Center,  the  wave  height  could  be  greater  than  five  feet  for  33 
percent  of  the  fall  and  winter  seasons  and  greater  than  eight  feet  for 
10  percent  of  the  same  seasons.  Under  these  conditions,  oil  spills 
during  transfer  operations  are  ”ery  likely  to  occur  and  environmental 
contamination  would  be  a serious  problem.  Before  developing  a new  | 

technology  to  counteract  potential  Oil-spill  problems,  efforts  must  be 
directed  toward  identifying  the  sea-state  threshold  beyond  which  con- 
tainment and  recovery  are  impractical  because  of  the  natural  dispersion 
of  oil  slicks.  The  determination  of  the  sea-state  threshold  for  oil 
spill  breakup,  however,  requires  an  understanding  of  the  turbulent 
motion  in  the  upper  ocean  and  its  effect  on  the  dispersion  of  oil  slicks. 

In  the  following,  let  us  review  our  understanding  of  this  phenomena. 

In  the  upper  ocean,  there  are  three  important  mechanisms  for  the 
generation  of  turbulence:  (i)  wave  breaking,  (li)  shear  or  current,  and 
(ill)  buoyancy-induced  motion.  The  Importance  of  the  first  mechanism, 
wave  breaking,  in  generating  turbulence  is  shown  by  the  measurements  of 
Grant,  et  al.  (1963),  i.e.,  the  turbulent  dissipation  rate  in  the  region 
near  the  surface  has  a larger  value  than  that  at  some  distance  below  the 
surface.  Phillips  (1966)  estimated  that  the  energy  flux  of  the  turbulence 
from  breaking  waves  may  be  comparable  with,  or  greater  than,  the  direct 
energy  flux  of  the  surface  current  associated  with  wind  stress.  The 
second  mechanism,  dynamic  instability  of  various  shearing  flows,  is 
mainly  caused  by  wind- Induced  currents.  Background  shearing  currents, 
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for  example,  caused  by  the  confluence  of  different  water  masses  are 
often  sources  of  appreciable  turbulence  in  surface  waters.  The  third 
mechanism,  buoyant  or  convective  Instability,  Is  important  in  regions 
where  substantial  cooling  or  surface  evaporation  appears  and  the  surface 
layer  tends  to  be  statically  unstable.  This  mechanism  results  In  free 
convection  characterized  by  an  irregular  cellular  motion. 

In  addition  to  the  three  mechanisms  of  turbulence  generation 
previously  described,  stratification  in  the  surface  water  requires  some 
consideration.  In  a stably  stratified  fluid,  stratification  can  inhibit 
the  development  of  turbulence  and  the  downward  diffusion  of  turbulent 
energy.  Another  aspect  of  the  stratification  effect  stems  from  the 
deepening  of  the  mixed  surface  layer,  and  thus,  the  thermocllne  below 
the  mixed  surface  layer.  A portion  of  the  energy  flux  from  the  surface 
is  used  to  deepen  the  thermocllne.  Hence,  the  dynamics  of  the  turbulence 
in  surface  waters  cannot  be  separated  from  those  of  the  mixed  layer. 

For  the  present  study,  a sea  in  a strongly  excited  state  is  of  concern 
and  thus,  the  stratification  In  surface  waters,  the  deepening  of  the 
mixed  layer,  and  the  bouyancy-lnduced  Instabilities,  are  expected  not  to 
be  Important  in  the  generation  of  turbulence.  Hence,  in  the  following 
paragraph  we  briefly  review  the  previous  works  concerning  breaking  waves 
and  wind-induced  currents. 

Breaking  waves  have  not  been  investigated  carefully  until  recently. 
Longue t-Higgina  (1969a,  1969b,  1973,  1974)  conducted  theoretical  studies 
on  the  mechanisms  of  wave  breaking.  Banner  and  Phillips  (1974)  investi- 
gated the  incipient  breaking  of  small-scale  waves.  Van  Dorn  and  Pazan 
(1975)  measured  the  velocity  in  breaking  waves  in  a wave  channel.  The 
wind- Induced  drift  currents  have  been  Investigated  by  many  researchers. 
Keulegan  (1951)  and  Plate  (1970)  measured  the  surface  current  in  wind- 
wave  tanks.  Baines  and  Knapp  (1965)  measured  the  drift  profiles  in  a 
shallow  channel.  Bye  (1965),  Shemdln  (1973)  and  Wu  (1975)  measured  the 
profiles  of  the  wind-induced  drift  current  In  wind-wave  tanks.  Turbulence 
measurements  in  the  water  near  the  ocean  surface  were  conducted  by 
Shonting  (1968),  by  Thorton,  et  al.  (1977)  and  by  Donelan  (1977). 


2 


In  order  to  understand  the  behavior  of  oil  dispersion  in  the  upper 
ocean,  it  is  necessary  to  know  the  turbulence  characteristics  (for 
example,  vertical  profiles  of  mean  velocity  and  turbulence  intensity). 

These  characteristics  are  not  yet  available  from  existing  literatures. 
Therefore,  the  first  objective  for  the  present  research  study  is  to 
conduct  laboratory  measurements  of  the  water  velocity  under  the  action 
of  wind  and  waves.  From  these  measurements,  one  could  obtain  an  empirical 
formula  so  as  to  apply  the  laboratory  results  to  the  field  situation. 

As  far  as  oil  dispersion  studies  are  concerned,  we  have  found  no 
systematic  work  on  oil  slick  breakup,  however,  some  works  on  oil  slick 
drift  are  available  and  they  are  summarized  briefly  as  follows. 

Schwartzberg  (1971)  conducted  laboratory  experiments  to  study  oil  spills 
under  the  action  of  wind.  He  measured  the  surface  drift  current  with  an 
oil  spill  present.  Pottinger  and  Reisbig  (1973)  conducted  laboratory 
measurements  of  the  drift  velocity  of  an  oil  lens.  Reisbig  et  al. 

(1973)  investigated  the  coupled  parallel  effects  of  wind  and  waves  on 
oil  spill  drift.  In  the  field,  Murray  (1975)  observed  the  effect  of 
wind  and  current  on  the  drift  of  large-scale  oil  slicks.  The  second 
objective  of  the  present  study  is  to  investigate  the  oil  dispersion.  We 
will  determine  the  vertical  penetration  and  the  size  distribution  of  the 
oil  droplets  when  the  oil  slick  is  broken. 

Our  approach  is  to  conduct  laboratory  experiments  in  a wind  tank. 
First,  we  investigate  the  water  turbulence  under  the  influence  of  wind. 
Second,  we  add  mechancially-generated  waves  in  the  tank  to  study  the 
effect  of  waves  on  the  water  turbulence.  And,  third,  we  select  experi- 
mental conditions  for  oil  dispersion  experiments. 

We  describe  the  experimental  setup  in  Section  2 and  the  experimental 
instruments  and  techniques  in  Section  3.  New  instruments,  such  as  the 
laser  Doppler  velocimeter  and  the  photodiode  wave  gauge  have  been  developed 
to  sense  remotely  the  velocity  and  the  vertical  displacement  of  the  water 
surface.  In  Section  4,  we  present  the  measurements  of  water  turbulence 
under  the  predominant  influence  of  wind.  In  Section  5,  the  effects  of 
waves  on  the  water  motion  are  presented.  Vertical  profiles  of  the  mean 
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velocity!  rns  velocity  fluctuations  and  other  statistical  quantities  are 
presented.  In  Section  6,  we  present  the  oil  concentration  measurements! 
Including  the  oil  concentration  profile  and  droplet  size  distribution. 
Conclusions  based  on  our  present  study  are  given  in  Section  7.  Experi- 
mental results  are  presented  in  tabular  and/or  figure  form.  Finally, 
recommendations  for  future  research  is  presented  in  Section  8. 
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2. 


EXPERIMENTAL  SETUP 


2.1  Wind-Wave  Tank 

The  wind-wave  tank  used  In  the  present  Investigation  is  composed  of 
a wave  tank  and  a wind  tunnel,  each  of  which  ie  9.1  m long,  1.2  a wide 
and  0.91  m deep.  Figures  la  and  lb  show  the  wind  tunnel  and  the  wave 
rank.  The  description  of  the  wind-wave  tank  has  been  presented  by 
Veenhuizen,  et  al.  (1975).  Since  Lne  wind-wave  tank  has  been  modified 
for  the  present  investigation,  a brief  description  is  given  as  follows. 

The  wind  tunqel  is  open-ended  with  a blower  located  in  the  upwind 
end.  The  wind  tunnel  is  suspended  from  above  over  the  wave  tank.  It 
has  an  adjustable  ceiling  to  allow  for  the  adjustment  of  the  pressure 
gradient  in  the  wind  tunnel.  The  maximum  speed  which  can  be  obtained  in 
this  wind  tunnel  is  10  m/sec. 

The  wave  tank  has  transparent  side  walls.  A carriage  was  mounted 
on  top  of  the  wave  tank  to  transport  instruments  and  to  support  probes 
within  the  wave  tank  and  the  wind  tunnel.  A wave  maker  was  installed  at 
the  upwind  end  of  the  water  tank.  The  wave  maker  is  a vertical  flap, 
made  of  20  mm  plexiglas  plate,  which  is  40  cm  high  and  spans  across  the 
wave  tame.  It  is  driven  by  an  electro-hydraulic  servo-control  and  is 
capable  of  generating  surface  waves  containing  single  or  multiple 
frequency  components. 

The  frequency  response  of  the  flap  displacement  is  essentially 
constant  over  the  frequency  range  of  1/2  to  10  Hz  with  virtually  no 
phase  shift.  The  maximum  stroke  of  the  flap  is  10  cm.  At  higher  fre- 
quencies, the  response  displacement  of  the  flap  diminishes  rapidly  and 
ia  limited  by  the  mechanical  response.  The  maximum  wave  height  is  about 
8 cm.  If  the  wave  height  is  over  8 cm,  appreciable  cross  waves  appear. 

A wave  absorber  was  located  at  the  downstream  end  of  the  tank  to 
absorb  surface  waves  generated  by  the  wind  or  the  wave  maker.  Basically, 
it  is  a rough,  porous,  sloping  beach  which  dissipates  the  surface-wave 
energy  by  causing  the  waves  to  break  as  they  run  up  on  the  beach  and  to 
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be  drained  of  energy  further  as  the  fluid  runs  back  down  the  beach 
through  many  small  openings  on  the  beach  surface.  The  wave  absorber 
(Figure  1)  is  constructed  of  a plastic  grid,  topped  with  5 cm  layers  of 
horse  hair,  and  covered  with  a layer  of  herring  fish  netting.  The  wave 
absorber  extends  across  the  width  of  the  tank  and  is  1.5m  long.  When 
the  absorber  is  at  about  a 10-degree  angle  to  the  water  surface  in  the 
tank,  the  reflection  coefficient,  defined  as  amplitude  reflected/amplitude 
incident,  is  measured  to  be  less  than  5 percent.  Thus,  the  wave  absorber 
works  effectively  in  damping  out  the  upcoming  waves. 

2 . 2 Reverse  Flow  Control 

When  wind  blows  over  the  water  surface,  a drift  layer  in  the  surface 
water  sets  up  due  to  the  shear  stress  exerted  on  the  v’ater  by  the  wind. 
Through  the  action  of  viscosity,  a water  boundary- layer  develops  which 
can  become  turbulent  when  the  Reynolds  number  is  sufficiently  large. 

The  boundary-layer  flow  entrains  the  fluid  from  below,  and  its  thickness 
grows  with  fetch.  To  satisfy  the  continuity  of  the  fluid  and  because  of 
the  finite  length  of  the  wave  tank,  a reverse  flow  results  in  the  lower 
part  of  the  tank.  The  drift  layer  and  the  reverse  flow,  together,  form 
a very  complicated  flow  field.  On  the  other  hand,  in  the  ocean  except 
the  near  shore  region,  such  a reverse  flow  is  negligible  due  to  the 
large  scale  of  the  flow  field.  Hence,  it  is  desirable  to  minimize  the 
reverse  flow  effects  in  the  laboratory  experiment. 

A false  bottom  15  cm  above  the  tank  floor  was  Installed  and  a tran- 
sition plate  guided  the  reverse  flow  close  to  the  wave  absorber  to  flow 
under  the  false  bottom.  In  the  upwind  end  of  the  false  bottom,  there 
was  a spacing  of  about  1 m between  the  end  wall  and  the  end  of  the  false 
bottom.  The  shape  and  dimensions  of  this  transition  plate  were  determined 
by  trial  and  error  to  maximize  the  flow  under  the  false  bottom  and 
minimize  the  reverse  flow,  if  any,  above  it.  For  the  case  when  wind 
blows  over  the  surface,  the  false  bottom  worked  quite  well  and  the  reverse 
flow  was  practically  nejligible  above  the  false  bottom.  When  mechanically- 
generated  waves  were  added,  the  false  bottom  did  not  seem  to  control  the 
reverse  flow  well,  so  it  was  removed.  The  measurement  here,  however. 
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showed  that  the  drift  layer,  for  the  case  of  mechanically-generated 

* 

waves  with  wind,  is  not  as  deep  as  for  the  case  when  wind  blows  over  the 
surface,  and  thus,  the  reverse  flow  control  may  not  be  critical. 

2.3  Oil  Feeding  System 

A flow  system  was  made  to  feed  a steady  and  uniform  layer  of  oil  at 

the  upwind  end  of  the  wave  tank,  and  to  remove  the  oil  continuously 

from  the  downwind  end  to  avoid  build-up  of  oil.  Figure  2 is  a schematic 

of  the  oil  feeding  system.  A Teel  rotary  gear  pump  manufactured  by 

Dayton  Electric  Manufacturing  Company,  model  ID  767,  capable  of  delivering 
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27  liter /min  with  no  head  and  25  llter/min  at  7 kg/cm  , is  used  to  move 
the  oil  from  a 230  litre  polyurethane  tank  to  a 20  litre  constant  head 
tank  that  is  about  1.8  m above  the  water  surface  in  the  wind-wave  tank. 

Oil  from  the  constant  head  tank  flows  through  a rotameter  to  a specially 
designed  oil  feeder.  The  feeder  consists  of  a 19  mm  aluminum  tube  with 
48  holes  ranging  in  sizes  from  3-  to  6 mm  and  pointing  upward.  Outside 
this  tube  is  a 30  mm  plastic  tube  with  48  holes,  5 mm  in  diameter,  25  mm 
apart  and  pointing  downward.  Outside  the  plastic  tube  is  a 32-  by 
32  mm  square-aluminum  channel  with  web  facing  downward.  A teflon  sheet 
with  one  end  taped  to  the  channel  guides  the  oil  flow  down  on  the  water 
surface.  The  other  end  of  the  teflon  sheet  floats  on  the  water  surface. 
This  oil  feeder  spans  the  tank  above  the  water  surface  and  below  a 
transition  plate  which  guides  the  air  flow  over  the  water  surface.  The 
oil  discharge  read  from  the  rotameter  was  calibrated  with  a container  of 
known  volume  «rd  a stop  watch.  It's  worth  mentioning  that  slots  were 
made  in  the  transition  plate  to  allow  air  to  enter  into  the  space  between 
the  transition  plate  and  the  wave  surface.  Thus,  the  teflon  sheet  was 
forced  to  be  flat  on  the  water  surface  and  to  follow  the  wave  motion. 

The  oil  fed  from  the  oil  feed  system  floats  on  the  water  surface  and 
drifts  downwind  by  the  current  Induced  by  wind  and  waves.  As  the  oil  slick 
reaches  the  downstream  end,  it  is  mixed  with  some  water  due  to  the 
turbulent  mixing.  The  oil-water  mixture  Is  then  removed  from  the  water 
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surface  using  a two-dimensional  overflow  drain  that  connects  to  a second 
230  liter  polyurethane  tank  which  serves  as  a settling  tank  to  allow  the 
oil  to  separate  from  the  water.  The  recovered  oil  cen  be  recycled  for 
the  next  experiment,  and  the  water  is  discarded.  This  semi-closed  system 
permits  a continuous  operation  for  about  15  minutes,  which  is  enough 
time  for  completion  of  an  experimental  run  of  oil  dispersion. 

2.4  Selection  of  Oil 

For  the  oil  dispersion  equipment,  it  is  usually  best  to  use  crude 
oils,  e.g.,  Kuwait  or  Libyan  crude.  We  did  not  use  crude  oils,  however, 
because  it  is  difficult  to  control  their  properties,  since  they  can 
vary  from  barrel  to  barrel.  Moreover,  the  purpose  of  the  present  study 
was  to  characterize  and  identify  dispersion  mechanisms  acting  on  an  oil 
slick  of  any  kind.  An  oil  was  selected  based  on  the  following  criteria. 
The  density  and  viscosity  of  the  oil  must  be  close  to  those  of  the  crude 
oils  mentioned.  The  oil  must  not  mix  easily  with  the  water  otherwise, 
the  oil  properties  could  change  significantly  as  the  dispersed  oil 
drifts  along  the  wind-wave  tank.  In  othar  words,  we  are  not  interested 
in  the  emulsification  process  occurring  in  the  present  experiments.  The 
oil  must  separate  easily  from  water  when  the  mixture  of  oil  and  water 
returns  to  the  oil  settling  tank  from  the  overflow  drain  located  in  the 
downwind  side  of  the  wave  absorber,  thus,  the  oil  can  be  recovered  for 
re-use.  The  oil  must  be  able  to  form  a thin  sheet  like  the  oil  slick  in 
the  ocean. 

To  examine  the  rate  of  recovery  and  separation  from  water  after 
mixing,  four  kinds  of  oils,  namely,  vegetable  oil,  Texaco  Almag  oil, 
Diesel  oil  No.  2 and  mineral  oil  were  tested  for  selecting  the  oil  to  be 
used  for  the  oil  dispersion  experiment.  The  oil  sample  was  put  in  a 
bottle  with  an  equal  amount  of  water  and  shaken  for  several  seconds.  A 
general  observation  as  to  the  rate  of  recovery  into  two  separate  fluids 
was  made  and  is  described  as  follows.  For  vegetable  oil,  the  rate  of 
recovery  is  slow.  Foam  appears  on  the  boundary  between  fluids,  and 
water  droplets  are  formed  and  are  suspended  in  the  oil.  For  Texaco 
Almag  oil  the  rate  of  recovery  is  also  slow,  and  water  droplets  sus- 
pended in  oil  are  observed.  Diesel  oil  No.  2 has  fast  recovery  and 
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clean  separation.  Mineral  oil  has  slow  recovery;  droplets  form  and 
suspend  in  both  fluids,  and  the  oil  clings  to  the  bottle. 

Spreading  on  a water  surface  also  was  tested  with  the  above  and 
additional  oils  in  selecting  the  best-suited  oil  for  the  oil  dispersion 
experiment.  A small  quantity  of  the  oil  sample  was  poured  on  a water 
surface.  The  resulting  oil  behavior  was  observed  and  is  described  as 
follows.  Diesel  Oil  No.  2 forms  a film  over  the  entire  surface.  Diesel 
Oil  No.  1 forms  a film,  but  not  as  readily  as  the  Diesel  Oil  No.  2. 

Mineral  Oil  forms  lenses,  and  a thick  film  adheres  to  the  glass.  Vege- 
table Oil,  Texaco  Almag  Oil  and  hydraulic  oil  No.  215  all  form  lenses. 

Based  on  the  tests  described  above,  we  selected  Diesel  oil  No. 

2 because  it  absorbs  water  least,  separates  from  water  quickly  and  can 
form  a thin  oil  film  over  the  water  surface.  Prior  to  each  experiment, 
the  physical  properties  of  the  oil,  such  as  specific  gravity  and  viscosity, 
are  measured.  Table  1 is  the  summary  of  the  ell  properties.  The 
physical  properties  did  not  vary  noticeably  during  the  period  of  experi- 
ments. For  example,  on  December  14,  1977,  we  measured  the  oil  properties 
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before  use,  i.e.,  viscosity  » 18.9  cmVsec  and  specific  gravity  * 0.84. 
After  use  during  the  same  day,  water  had  mixed  with  oil,  and  both  viscosity 
and  specific  gravity  increased  accordingly.  New  oil  of  about  a quarter 
volume  of  the  used  oil  was  then  added  to  the  uBed  oil  to  ensure  that  the 
oil  properties  remained  practically  unchanged.  For  example,  between 
December  14,  1977  when  new  oil  was  first  used  and  January  18,  1978, 
viscosity  and  specific  gravity  increased  by  about  7-  and  10  percent 
respectively.  Whenever  we  found  that  the  specific  gravity  and  viscosity 
changed  significantly,  the  recycled  oil  was  drained  and  replaced  by 
fresh  oil.  As  noted  in  Table  1,  new  oil  was  supplied  for  future  experi- 
ments on  January  19,  1978.  The  interfacial  tension  at  20°  C between 
air  and  water  and  between  air  and  oil  are  about  70-  and  30  dyne/cm, 
respectively. 
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EXPERIMENTAL  INSTRUMENTS  AND  TECHNIQUES 


3. 1 Hot-Wire  and  Hot-Film  Probee 

For  velocity  measurements  In  the  air,  we  used  two  hot-wire  probes 
manufactured  by  Thermo  Systems  Inc.  (TSI  1210-T1.5)  operated  with  two 
constant-temperature  hot-wire  anemometers  (TSI  Model  1054B).  A typical 
calibration  curve  of  the  hot-wire  probe  is  plotted  in  Figure  3.  A 
Pitot-static  tube  with  a Gilmont  micromanometer  was  used  for  calibration. 

Twenty  channels  of  hot- film  anemometer  (TSI  Models  1051-10  and 
1053B)  were  used  with  ten  cross-film  probes  (TSI  Model  1284)  to  make 
simultaneous  measurements  of  the  instantaneous  longitudinal  and  vertical 
velocity  components  in  the  water.  The  hot- film  probes  were  mounted  on  a 
probe  strut  which  in  turn  was  clamped  on  the  instrument  carriage.  The 
calibration  was  made  by  towing  the  hot-film  probes  along  the  wave  tank. 
The  towing  speed  was  registered  with  the  rotational  speed  of  the  motor 
which  drove  the  towing  system.  A slot  wheel  with  a light-emitting-diode 
and  a photocell  was  mounted  on  the  motor  6haft.  As  the  motor  turns,  the 
photocell  registers  a signal  which  corresponds  to  the  speed  at  which  the 
slot  wheel  and  the  motor  rotate.  Figures  4(a)  and  (b)  show  two  cali- 
brations of  a hot-film  probe  for  high  and  low  speed  ranges. 

3.2  Laser  Doppler  Veloclmeter 

The  hot-film  probe  is  an  accurate  and  sensitive  instrument  for 
measuring  velocity;  however,  if  there  exists  a reverse  flow,  the  hot- 
film  probe  is  not  useful  because  it  senses  speed  but  not  the  flow 
direction.  Hence,  we  cannot  use  the  hot-film  probe  to  measure  velocity 
in  large  amplitude  waves,  where  large  reverse  flow  results  from  the 
oribital  motion  of  the  waves.  To  overcome  this  difficulty,  we  use  a 
laser  Doppler  veloclmeter  (LDV)  developed  at  the  Flow  Research  Company. 
The  LDV  was  designed  to  operate  in  the  dual-beam  backscatter  mode,  and 
it  has  the  frequency  shift  capability  to  measure  velocity  even  when  the 
reverse  flow  is  present.  Some  specifics  which  are  relevant  to  the 
present  experiments  are  described  briefly  as  follows.  Detailed  descrip- 
tions of  the  design  are  given  by  Bachalo  and  Lin  (1977). 
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The  LDV  consists  of  a 4-watt  Argon-ion  laser  (Spectra  Physics  Model 
164-05),  a set  of  transmitting  and  receiving  optics,  two  acousto-optic 
modulators  or  Bragg  cells  (Isomet  Corporation),  a photomultiplier  (RCA 
7767/1C  with  Sll  response),  an  amplifier  (Advantek  SD7-0614B) , a counter- 
type signal  processor  (Macrodyne  Model  2096)  and  a bandpass  filter 
(Krohn-Hite  Model  3103).  The  laser,  the  optics  and  the  photomultiplier 
were  all  mounted  on  a 1.2-  by  0.91  m optical  table  which  was  made  of 
0.64  cm  aluminum  plates  sandwiched  with  honeycomb.  The  top  surface  of 
the  optical  table  was  drilled  with  0.64  cm  threaded  holes  into  which 
optic  stands  can  be  mounted.  The  optical  cable  was  mounted  on  a steel 
base,  which  was  equipped  with  an  oil  jack  for  adjustment  of  vertical 
displacement  up  to  48  cm,  The  surface  also  could  be  tilted  in  one 
direction  about  + 30  deg.  from  its  leveled  position. 

Figures  5 (a)  and  (b)  are  a photograph  and  a schematic  of  the  LDV.  The 


laser  bean,  projected  from  the  laser  housing  is  reflected  twice  by  two 
45  deg.  mirrors,  and  then  9plit  into  two  beams  with  a splitter  cube. 

Both  beams  were  frequency  shifted  with  two  Bragg  cells  and  were  then 
focused  with  an  achromatic  lens  at  the  measuring  point.  Due  to  the 
interference  of  the  two  spilt  light  beams,  a fringe  pattern  forms  at  the 
focal  point,  and  the  fringe  pattern  moves  at  a speed  corresponding  to 
the  frequency  shifted  by  the  Bragg  cells.  The  speed  at  which  the  fringe 
frame  moves  is  calculated  by  x M , where  d^  is  the  fringe 
spacing  and  Af  is  net  frequency  shift.  The  fringe  spacing  d^  can  be 
calculated  by 

df  - a (3-1) 

2 sin  (|) 

where  X is  the  wavelength  of  the  laser  and  0 is  the  beam  crossing 
angle  at  the  focal  point. 

The  water  was  seeded  with  Pliolite  plastic  particles  for  scattering 
purposes.  The  particles,  manufactured  by  Goodyear  Rubber  Company,  were 
ground,  sieved  and  filtered  to  remove  large  ones.  Figure  6 shows  a photo- 
graph of  the  particles  viewed  with  an  electronic  microscope.  The  average 


■•4 

1 


<4 


i 


1 1 


particle  size  is  less  than  40  y.  The  particles  had  a specific  gravity 
of  about  1.03.  Because  of  their  small  size,  the  particles  can  follow 
the  flow  motion  very  well. 

A receiving  lens  (0.5  m focal  length)  imaged  the  focal  point  onto  a 
pnotomultiplier  through  a pinhole  0.05  cm  in  diameter.  The  scattering 
light  from  a seeded  particle  at  the  focal  point  was  therefore  collected 
onto  the  photomultiplier.  The  photomultiplier  output,  which  was  propor- 
tional to  the  scattered  light  intensity,  was  amplified,  band-pass 
filtered  and  fed  into  the  Macrodyne  signal  processor. 

The  Macrodyne  processor  can  process  and  validate  the  frequency 
information  in  the  signals,  and,  then,  provide  voltage  output  in  both 
analog  and  digital  forms  proportional  to  the  Doppler  shift  period  T ^ . 
Figure  7 is  a block  diagram  showing  the  various  stages  of  data  recording 
and  processing  of  the  LDV.  The  particle  velocity  Up  can  thus  be 
calculated  from  the  difference  between  the  Doppler  shift  frequency  and 
the  frequency  shift  introduced  by  the  Bragg  cells, 

Up  - df(fd  - Af)  (3-2) 

where  f^  » 1/t^  is  the  Doppler  shift  frequency.  The  maximuT  negative 
velocity  that  can  be  resolved  by  the  LDV  is  therefore  U = -dp  Af  . 

The  setup  can  be  slightly  modified  for  measurements  of  the  instantaneous 
vertical  velocity  by  rotating  the  bean  splitter  cube  and  the  two  Bragg 
cells  90  deg.  Figure  8 shows  a calibration  curve  of  the  LDV  system 
measured  in  an  air  Jet.  For  the  measurements  in  the  air,  we  used  water 
droplets  as  seeding  particles.  The  LDV  measurement  is  in  good  agreement 
with  the  Pitot  tube  result.  The  calibration  of  the  LDV  in  water  is 
described  in  Appendix  B. 

3.3  Capacitance  Probe 

For  wave-height  measurment,  we  used  a capacitance  gauge  manufactured 
by  Flow  Research  Company  (Model  1204).  The  capacitance  gauge  had  a bare 
stainless  steel  tube  (1.6  mm  o.d.)  or  a teflon  coated  wire  as  the  sensor. 
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Calibration  was  made,  before  or  after  each  run,  by  displacing  the  probe 
at  several  predetermined  vertical  increments.  Figure  9 is  a typical 
calibration  curve  of  the  capacitance  probe.  The  calibration  points  were 
fitted  with  a second  degree  polynomial. 

3.4  Photodiode  Wave  Gauge 

In  case  an  oil  slick  is  present,  the  capacitance  probe  falls  to 
measure  surface  wave  height  because  of  oil  contamination.  To  avoid 
contamination,  we  developed  a photodiode  wave  gauge  - a remote  sensor 
for  measuring  surface  wave  height.  The  photodiode  wave  gauge  includes  a 
Reticon  camera  and  a controller  manufactured  by  Reticon  Corp.  (Model  LC 
600  V and  RS  605)  and  a 4-watt-Argon-ion  laser  (Spectra  Physics,  Model 
164-05).  Figure  10  shows  a schematic  of  the  photodiode  wave  gauge, 
and  Figure  11  is  a photograph  of  the  experimental  setup  for  remotely 
sensing  surface  waves.  The  laser  beam  serves  as  a light  source  to 
illuminate  the  interface  between  the  air  and  the  water  or  between  the 
air  and  the  oil.  Fluorescent  dye  (Fluorescein  disodium  salt)  was  pre- 
mixed uniformly  in  the  water.  The  concentration  was  about  1 ppm.  By 
vertically  projecting  the  laser  beam  from  the  air  to  the  water  or  to  the 
oil  slick,  a bright  fluorescent  beam  is  produced  in  the  water  and  a 
sharp  dark-to-brlght  interface  appeared  on  the  water  surface  or  the 
surface  of  the  oil  slick.  This  dark-to-hright  interface  can  be  detected 
by  a vertical  linear  array  of  photodiodes  mounted  inside  the  Reticon 
camera. 

The  photodiode  array  has  256  elements  spaced  25.4  ym  center  to 
center.  The  air-water  interface  marked  with  the  fluorescent  dye  was 
focused  on  the  photodiode  array  by  optics  including  lenses  and  extension 
tubes.  The  f leld-of-view  can  vary  from  a fraction  of  a centimeter  to  a 
few  meters.  The  spatial  resolution  in  measuring  the  vertical  displacement 
of  the  air-water  interface  can  be  adjusted  by  varying  the  size  of  the 
field-of-view.  For  example,  the  spatial  resolution  is  0.04  cm  when  the 
field-of-view  is  10  cm. 

The  position  of  the  image  of  the  air-water  interface  cr  air-oil 
Interface  was  determined  electronically  by  a scanning  control.  The 
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scanning  rate  ranges  from  0.4-  to  40  ys,  and  thus,  it  ensures  the  fast 
response  of  the  photodiode  wave  gauge.  Figure  12a  shows  a schematic  of 
the  laser  beam  and  the  corresponding  light  intensity  profile  along  the 
beam  in  the  vicinity  of  the  air-water  interface.  Figure  12b  shows  a 
complete  scan  of  the  VIDEO  output  of  photodiodes.  The  scan  started  from 
the  first  diode  which  was  focused  on  the  laser  beam  above  the  air-water 
interface  and  ended  at  the  last  diode  which  was  focused  on  the  laser 
beam  below  the  interface.  The  vertical  lines  in  Figure  l 2b  represent 
the  analog  pulses  of  the  individual  diodes,  and,  at  the  interface  there 
is  a large  jump  in  voltage  output.  When  surface  waves  are  present,  the 
air-water  interface  oscillated  vertically  and  the  pattern  of  the  analog 
pulses  shown  in  Figure  12b  moves  accordingly.  By  detecting  the  position 
of  the  sharp  slope  of  the  analog  pulse  one  could  therefore  register  the 
vertical  displacement  of  the  water  or  oil  surface.  The  registration  was 
done  by  comparing  a THRESHOLD  output  with  the  VIDEO  output.  When  the 
voltage  output  of  an  analog  pulse  is  Jess  than  the  threshold  level,  the 
processor  registers  a logical  "zero".  When  the  voltage  output  exceeded 
the  threshold  level  within  one  ecan,  the  processor  will  register  logical 
"ones"  for  the  rest  of  the  scan  (Figure  12c).  The  registered  logical 
signal  was  then  electronically  counted  to  determine  the  position  of  the 
"nth"  diode  where  the  "zero"-to-"one"  occurred  and  where  the  air-water 
interface  is  imaged. 

Calibration  of  the  photodiode  gauge  was  made  by  displacing  it  at 
several  vertical  positions  with  predetermined  increments.  Die  diode 
numbers  are  plotted  versus  the  vertical  displacement  in  Figure  13.  The 
calibration  points  are  fitted  with  a straight  line. 

The  performance  cf  the  photodiode  wave  gauge  was  dynamically  cali- 
brated with  a capacitance  probe  positioned  side  by  side  in  close  proximity 
in  a mechanically-generated  wave.  The  wave  frequency  was  selected  to  be 
1.6  Hz  for  which  the  meniscus  effects  cn  the  performance  of  the  capaci- 
tance probe  are  not  important.  In  Figures  14  and  15,  we  present  the 
wave  profiles  of  mechanically-generated  waves  measured  with  the  photodiode 
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wave  gauge  and  Che  capacitance  probe.  The  two  wave  profiles  agree  quite 
well  except  in  the  regions  near  wave  crests  and  troughs.  Apparently, 
near  wave  crests  and  troughs,  the  surface  tension  effect  causes  the 
deviation  between  the  wave  profiles  measured  with  the  capacitance  probe 
and  the  photodiode  wave  gauge.  Note  that  Sturm  and  Sorrell  (1973) 
conducted  an  experiment  to  determine  the  spatial  and  frequency  response 
of  a capacitance  gauge  to  surface  waves.  Because  of  meniscus  effects 
due  to  surface  tension,  the  capacitance  probe  has  a spatial  resolution 
of  about  6 mm  and  a frequency  response  of  6 Hr. 

3 . 5 Data  Acquisition  and  Process 

The  signals  measured  with  instruments  described  in  previous  sections 
are  transmitted  to  a minicomputer  system  which  stores  and  processes 
the  data.  The  system  consists  of  a Nova-800  minicomputer  with  a 32-K 
core,  suported  by  peripheral  equipment.  This  peripheral  equipment 
includes:  (1)  two  magnetic  disk  drives,  an  Iomec  and  a Diablo,  each 
having  a 2.5  M words  capacity,  (2)  a Wang  9-track  800-BPI,  magnetic- 
tape  drive,  (3)  a Houston  Instruments  DPI  incremental  plotter,  and  f4)  a 
Ver6atec  model  1100A  electrostatic  printer-plotter.  In  addition,  the 
computer  is  equipped  with  floating-point,  integer-multiply-di vide , 
digltal-I/O,  and  digital-analog-conversion  hardware. 

Analog  data  from  the  wind-wave  tank  is  sent  through  signal  condi- 
tioners which  adjust  the  voltage  ranges  by  subtracting  a d.c.  level, 
applying  an  amplification  factor,  and  then  low-pass  filtering  the  signals 
before  analog-to-digltal  conversion  (A/D)  for  the  computer.  Signal 
conditioning  io  provided  for  48  channels  of  analog  data.  The  A/D 
converter,  Analogic  Corporation  Series  AN5800,  has  a capacity  of  48  data 
channels  at  13-bit  resolution,  and  an  overall  sampling  rate  of  44,000 
saroples/see . The  data  are  multiplexed  and  sent  to  the  computer  interface 
by  the  A/D  converter.  The  computer  records  the  digital  data  on  a magnetic 
disk  on-line  during  the  progress  of  the  experiments  for  playing  back 
and  processing  at  a later  time.  Data  stored  on  the  disk  may  be  transferred 
to  magnetic  tape  for  long-term  storage,  or  it  may  be  called  Into  the  com- 
puter core  for  processing. 
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4. 


EFFECT  OF  WIND  ON  WATER  MOTION 


Measurements  of  velocity  In  the  water  were  conducted  for  the  case 
when  the  wind  dominated  the  motion.  As  waves  were  generated  by  wind, 
the  wave  effect  on  the  water  motion  was  also  investigated.  We  first 
measured  the  mean  velocity  and  the  rms  longitudinal  velocity  in  the  air 
boundary-layer  over  the  water  surface.  Based  on  these  measurements,  we 
could  establish  the  wind  conditions  required  for  the  present  investi- 
gation. The  water  surface  displacement  was  measured  to  obtain  the  wave 
characteristics  for  the  experimental  conditions  which  were  of  interest 
here.  The  drift  velocity  on  the  water  surface  was  measured  to  determine 
the  effect  of  wind  speed  on  the  surface  drift,  velocity.  The  motion  in 
the  water  boundary-layer  was  visualized  with  dye,  and  the  boundary-layer 
thickness  was  determined  from  the  visualization  results.  The  effect  of 
wind  speed  and  fetch  on  the  evolution  of  the  boundary-layer  thickness 
was  established  from  the  visualization  results.  A vertical  array  of 
hot-film  probes  was  positioned  along  the  centerline  of  the  wind^wave 
tank  at  several  fetches  to  measure  the  velocity  field  in  the  water 
boundary-layer.  Vertical  profiles  of  the  mean  velocity,  the  rms  longi- 
tudinal and  vertical  velocity  fluctuations,  the  Reynolds  stress  and 
the  turbulence  dissipation  rates  were  obtained.  Velocity  spectra  were 
calculated  from  the  data  of  velocity  fluctuations  to  estimate  the 
contribution  of  the  wave  motion  on  the  velocity  field.  From  the  measure* 
ments  to  be  described  in  the  following,  we  could  then  determine  the 
effect  of  wind  on  the  motion  in  the  water  boundary-layer  under  a mobile 
water  surface. 

A definition  sketch  for  the  air  and  water  boundary-layers  to  be 

investigated  here  is  drawn  in  Figure  16.  The  air  boundary- layer  ia 

characterized  by  a thickness  6 , and  the  water  boundary-layer  by  a 

a 

thickness  6 . The  fetch,  x , is  measured  from  the  downstream  end  of 

w 

the  transition  plate  which  guides  the  air  flow  over  the  water  surface, 
and  z is  positive  upwards  measuring  from  the  undisturbed  water  surface. 
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The  dominant  wavelength  and  height  are  denoted  by  X and  H . The 
freestream  velocity  la  denoted  by  . 

4. 1 Air  Boundary-Layer 

Figure  17  shows  the  mean  velocity  and  turbulence  intensity  profiles 
measured  In  the  air  boundary-layer  at  x * 5 m and ’ - JO  m/sec.  The 
profiles  are  similar  to  the  ones  measured  by  others  in  wind-wave  tanks. 

The  rm8  of  longitudinal  velocity  fluctuations  has  a maximum  of  about 
12  percent  of  the  freestream  wind  speed.  The  boundary- layer  thickness 
Is  about  18  cm.  A logarithmic  profile  is  observed  when  plotting  the 
mean  velocity  profile  on  a semi-log  plot,  and  a solid  line,  representing 

^/u*a  " 5-75  lo8  i/z0  * 

is  drawn  on  Figure  18  to  approximate  the  logarithmic  profile.  The 
friction  velocity  estimated  from  best  fit  is  u*fl  • 58  cm/aec  and  the 
aerodynamic  roughness  16  zq  - 0.04  cm.  The  maximum  rms  of  longitudinal 
velocity  fluctuation  1b  about  twice  that  of  the  air  friction  velocity. 

These  results  show  that  the  air  turbulent  boundary- layer  over  wind-waves 
is  very  similar  to  the  turbulent  boundary- layer  over  a rough  flat  plate, 
and  they  are  consistent  with  previous  measurements  by  many  investigators. 

For  a turbulent  air  boundary- layer  over  wind-wavas,  Lin  and 
Veenhulzen  (1975)  collected  experimental  results  of  many  Inver tigators . 

They  correlated  the  growth  of  the  air  boundary-lcyer  thickness  with 

the  freestream  velocity  and  the  letch  by  the  expression 

g j 

- (0.025+0.003)  (*f)4/5,  for  0.2  < **  < 17  , (4-1) 

U U U 

OO  CO  oo 

where  g is  the  gravitational  accelaratlon.  They  also  found  that 
the  air  friction  velocity  u*  can  ba  axpressad  in  the  following  form 
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(4-2) 


r~—  - (0.055+0.01)  (^  )"1/10  , for  0.3  < M < 200  . 

u u± 

ao  cd 

The  present  data  conform  closely  to  these  empirical  formulas. 

4 . 2 Wave  Characteristics 

A typical  output  of  a capacitance  probe  is  shown  in  Figure  19. 

The  freestream  wind  speed  was  10  m/sec  and  the  capacitance  probe  was 

placed  at  5 m fetch.  The  Instantaneous  water  surface  displacement  <1 

is  measured  from  the  undisturbed  water  surface  level.  In  this  figure,  a 

dominant  frequency  is  evident,  and  wave  packets  of  smaller  frequency  are 

shown  to  be  consistent  with  other  wind-wave  observations.  The  frequency 

spectrum  of  the  water  surface  displacement  (wave  energy  spectrum), 

part  of  which  is  shown  in  Figure  19,  is  presented  in  Figure  20.  The 

dominant  frequency  is  3.5  Hz,  and  the  spectrum  for  4<f<10  hz  follows 
2 -5  -4 

S(f)  « a g f , where  a ■ 5 x 10  , for  this  particular  run.  The 

value  of  a is  larger  than  the  value  a - 7.5  x 10  ^ observed  in  the 
ocean  (Phillips,  1966).  Others  have  found  that  larger  values  of  a 
are  observed  in  the  laboratory  than  in  the  ocean  also.  The  spectrum  for 
frequencies  above  10  Hz  should  be  viewed  with  caution  since  the  capacitance 
probe  is  limited  in  its  frequency  response  due  to  the  meniscus  effects. 

The  variation  of  the  root-  mean-square  value  of  the  Water  surface 
displacements  with  wind  speed  Is  shown  in  Figure  21  for  two  fetches, 
x ■ 5-  and  6 m.  As  expected , v n increases  with  both  fetch  and  wind 
speed.  The  effect  of  wind  speed  variation  on  the  dominant  wave  frequency 
is  shown  in  Figure  22  for  two  fetches,  x ■ 5-  and  6 m.  Here,  the 
dominant  wave  frequency  f^  is  defined  as  the  frequency  at  which  the 
frequency  spectrum  peaks.  The  dominant  wave  frequency  decreases  with 
both  fetch  and  wind  speed.  The  effect  of  fetch  on  the  dominant  wave 
frequency  is  shown  in  Figure  23  for  wind  speed  Um  * 10  m/sec.  End  wall 
effects  in  the  wave  tank  possibly  contribute  to  the  leveling  off  of  the 
curve  shown  for  fetch  x > 5 m. 

The  dominant  wave  phase  speed  could  be  determined  by  correlating 
the  signals  from  two  capacitance  probes  separated  by  a distance  Ax  in  the 
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wind  direction.  For  a simple  train  of  wavas,  the  phase  speed  can  be 
Ax 

measured  as  Cp  " » where  A*  ie  the  time  lag  at  which  the  cross- 

correlation  has  its  first  maximum.  For  a random  train  of  waves,  as  is 
the  case  with  wind-waves,  a more  accurate  determination  of  At  for  each 
wave  component  requires  computing  the  Fourier  transform  of  the  croas- 
correlation.  The  complex  Fourier  transform  consists  of  a real  part  (co- 
spectrum)  and  an  imaginary  part  (quad-spcctrum) . At  the  dominant  frequency 

At  is  related  to  the  phase  angle  (■  tan  ^ c.°.  by 

quad-spectrum  ’ 

Phs36  dnslc 

At  “ — - j; -r6 — . Figure  24  shows  the  effect  of  wind  speed  variation  on 


the  wave  phase  speed  Cp  for  two  fetches. 


x » 5 and  6 m.  The  phase 
speed  Increases  with  increasing  wind  speed.  Figure  25  shows  fetch 


effects  on  the  wave  phase  speed  C 


The  end  wall  effects  possibly 


contribute  to  the  leveling  off  of  the  curve. 

ie  dominant  wave  frequency  f a 
one  can  easily  compute  the  dominant  wavelength 


From  the  dominant  wave  frequency  f and  the  dominant  wave  phase 


speed  C_ 

Cp  1 

X - t~  . Figure  26  shows  the  effect  of  wind  speed  variation  on  X for 
ro 

two  fetches,  x * 5 and  6m;  X increases  with  both  fetch  and  wind 
speed.  The  dominant  wave  slope  \y n /X  is  shown  in  Figure  27  and  it 
increases  with  fetch  and  wind  speed.  Figure  28  shows  the  effects  of 
fetch  on  the  dominant  wavelength  for  wind  speed  = 10  m/sec. 


4 . 3 Surface  Drift  Velocity 

The  drift  velocity  on  the  water  surface  was  measured  using  paper 
disks  which  had  been  coated  lightly  with  wax  so  they  would  float.  The 
diameter  of  the  disk  was  about  7 mm,  and  the  trajectory  of  the  disc 
floating  on  the  water  surface  was  recorded  with  a 16-mm  movie  camera. 

The  time  required  to  travel  a known  distance  was  registered  with  a 
digital  clock  with  0. 1 second  resolution.  For  each  run,  the  measurements 
were  repeated  several  times,  and  the  mean  and  the  standard  deviation  of 
the  measurements  were  determined.  Table  2 lists  the  measured  surface 
drift  velocity  in  a tabulated  form.  Also  shown  are  the  number  of  floats 
observed  for  each  run  as  well  as  the  standard  deviation.  Figure  29 
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shows  the  variation  of  surface  drift  velocity  Ufl  with  fetch  for  wind 
speeds  * 4,  7 and  10  m/sec.  The  surface  drift  velocity  Ug  has 
an  average  value  of  about  3.2  percent  of  the  wind  speed. 

4.4  Water  Boundary-Layer 

The  water  boundary-layer  that  develops  under  the  action  of  wind  was 
visualized  with  the  dye  technique.  A flourescent  dye  was  injected  at  the 
beginning  of  the  test  section  using  an  L-shape  tube  2 mm  in  diameter. 

A sheet  of  laser  light  was  used  to  illuminate  the  fluorescent  dye  in  a 
vertical  plane  parallel  to  the  wind  direction,  and  a movie  film  recorded 
the  motion  in  the  boundary  layer.  Figures  30  (a)  through  (f)  show  the 
photographs  of  the  dyed  boundary  layer  for  three  wind  speeds  ' U ■ 4-, 

7-  and  10  m/3ec  and  two  fetches,  x ■ 1.7-  and  1.2  n.  In  the  photographs, 
a bright  horizontal  line  appears  which  represents  the  water  surface. 

The  bright  line  is  produced  when  the  sheet  of  laser  light  intersects  the 
water  surface.  For  ■ 4 m/sec,  the  boundary  layer  appears  to  be 
laminar  at  x ■ 1.7  m,  and  it  becomes  turbulent  at  * * 3 a as  indicated 
by  the  substantial  increase  in  thickness.  For  - 7-  and  10  m/sec, 
the  flow  is  turbulent  at  both  fetches  as  indicated  by  irregularity  of 
the  dye  boundary  and  by  the  eddy  structure  in  the  boundary  layer.  From 
the  dye  visualization  results,  we  determined  the  mean  thickness  and  its 
standard  deviation  of  the  water  boundary-layer.  Table  3 contains  the 
results,  and  Figure  31  shows  the  dependence  of  the  water  boundary-layer 
thickness  on  fetch  and  wind  speed.  The  boundary-layer  thickness  Is 

non-dimenslonalized  using  the  wind  speed  and  the  acceleration  of 

gravity  g . The  solid  line  on  the  figure  represents  the  empirical 
correlation 

- 0.025  (-^y)^5  for  0.2  < *y  < 3 . (4-3) 

u u u 

OQ  OO  00 

The  coefficient  expressed  in  equation  4-3  is  coincidentally  equal  to 
that  computed  by  Lin  and  Veenhulzen  (1975)  for  air  boundary-layers  in 
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wind-wave  tank  experiments.  The  present  data  at  small  fetches  for  wind 
speed  Uw  • 4 m/sec  deviate  considerably  from  this  formula.  We  observed 
visually  that  at  these  fetches  the  water  boundary- layer  at  wind  speed 
•*  4 m/sec  did  not  resemble  the  fully-developed  turbulent  boundary- 
layers  that  were  observed  at  the  same  fetches  for  higher  wind  speeds. 


4 . 5 Mean  Velocity 

Mean  velocity  profiles  in  the  water  boundary-layer  were  obtained  by 
averaging  the  signals  from  individual  cross-element  hot-film  probes  over 
a period  of  40  seconds.  A typical  mean  velocity  profile  is  shown  in 
Figure  32  for  wind  speed  ■ 10  m/sec  and  fetch,  x - 5 m.  The 
figure  shows  three  different  runs  to  check  repeatability.  At  the  water 
surface,  the  mean  drift  velocity  determined  from  the  floating  disk 
measurement  is  plotted,  and  it  is  denoted  by  a cross  symbol.  Near  the 
water  surface,  the  data  show  a sharp  velocity  gradient.  The  drift 
layer  determined  from  the  mean  velocity  profile  has  a thickness  of  about 
25  cm.  Figure  33  is  a replot  of  the  mean  velocity  defect  profile  on  a 
semi-log  paper.  The  mean  velocity  is  shown  relative  to  the  surface 
drift  velocity.  One  could  fit  a straight  line  to  the  part  of  the  profile 
close  to  the  water  surface  by 

U - U 

-2 - 5.75  log  (.—■)  . 

u . z 


From  the  plot,  the  friction  velocity  u*  is  determined  to  be  1.81  cm/sec, 
and  the  roughness  height  zq  ■ 0.029  cm.  Thus,  the  present  result  shows 
the  existence  of  the  logarithmic  profile  of  the  mean  velocity  in  the 
region  near  the  water  surface,  and  it  is  consistent  with  previous  meaeure- 
t/.f>rts  in  wind  waves  by  Shamdin  (1973)  and  Wu  (1975).  The  measurement  of 
the  mean  velocity  defect  shows  that  the  water  boundary- layer  under  the 
water  surface  resembles  the  turbulent  air  boundary-layer  over  the  water 
surface,  and  the  water  boundary- layer  is  produced  as  the  result  of  the 
wind  stress  exerted  on  the  water  surface. 


Table  4 shows  the  friction  velocity  as  determined  from  similar 

profiles  for  different  wind  speeds  and  fetches.  Also  shown  is  the 

corresponding  friction  velocity  for  the  air  boundary-layer  as  computed 

from  the  empirical  formula  expressed  in  Equation  (4-2).  The  ratio  of 

2 2 

the  water  shear  stress  (pu^)w  to  the  air  shear  stress  (pu^)a  at  the 

air-water  interface  Indicates  the  percentage  of  momentum  that  is  extracted 

from  the  air  boundary-layer  by  the  water  boundary-layer.  This  ratio 

increases  as  the  fetch  and  the  wind  speed  increases.  At  U ■ 10  m/sec, 

and  x * 6 m , the  ratio  has  a value  of  0.9,  or  equivalently,  the  ratio 

of  the  wave  drag  coefficient  to  the  total  wind  drag  coefficient  is  0.1. 

Based  on  the  saturated  wave  energy  spectrum,  Phillips  (1966)  estimated 

the  ratio  of  the  wave  drag  coefficient  to  the  total  wind  drag  coefficient 

c /c,  to  be  about  0.1;  also  Lin  and  Veenhuizen  (1975)  obtained  the 
v a 

ratio  c /c,  for  wind  waves  by  correlating  many  experimental  data 
v a 

obtained  in  wind-wave  tanks,  i.e.. 


cw  m 3 pw  g n2 
Cd  8 Pa  U2  0 


(4-4) 


where  p^  and  p^  are  water  density  and  air  density  and  0 is  the 

momentum  thickness  of  the  air  boundary-layer.  Their  computed  ratio 

was  c / c,  * 0.1  + 0.03,  for  0.3  < -rr  < 17  . 
w d - - u2  * 

00 

Figure  34  shows  mean  velocity  profiles  in  the  water  boundary-layer 
at  fetch  x ■ 5 m,  for  three  different  wind  speeds,  U ■ 4-,  7-  and 
10  m/sec.  The  average  velocity  gradient  and  the  thickness  of  the 
boundary-layer  both  increase  with  wind  speed. 


4.6  RMS  Velocity  Fluctuations 

Figure  35  shows  the  vertical  profiles  of  the  rms  value  of  longi- 
tudinal velocity  fluctuations  in  the  water  boundary-layer  for  fetch 
x " 5 m and  wind  speed  U • 10  m/sec  for  three  repeated  runs.  The 
rms  values  monotonically  decrease  with  depth  below  the  water  surface. 

Its  maximum  value  occurs  close  to  the  water  surface  and  is  about  20  percent 
of  the  mean  surface  drift  velocity  and  about  3.5  times  the  frictional 
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velocity  u^  in  .water.  In  the  turbulent  air  boundary-layer  over  the 
water  surface,  previously  measured  by  Young  et  al.  (1973)  end  Yu  and  Lin 
(1973),  the  maximum  rma  longitudinal  velocity  fluctuations  is  about  IS 
percent  of  the  air  freestream  velocity  and  is  about  three  times  the  air 
frlccion  velocity.  Hence,  we  have  more  evidence  supporting  that  the 
water  boundary-layer  under  the  water  surface  resembles  the  air 
boundary-layer  over  the  water  surface.  It  should  be  noted,  however,  that 
contains  in  part  the  contribution  from  the  orbital  motion  of  the 
waves.  It  will  be  shown  in  Section  6.7  that  this  contribution  is  by  no 
means  negligible. 

Figure  36  shows  the  vertical  profiles  of  the  longitudinal  velocity 

fluctuations  measured  at  fetch  x ■ 5 m and  for  three  different  wind 

speeds  U ■ 4-,  7-  and  10  m/sec.  The  maximum  rms  value  is  about  20  percent 

of  the  mean  surface  drift  velocity  for  wind  speeds  7-  and  10  m/sec,  and 

10  percent  for  wind  speed  4 m/sec.  These  results  show  that  the  water 

boundary-layer  for  Uw  - 4 m/sec  has  not  been  fully  developed  when 

comparing  with  those  for  » 7-  and  10  m/sec. 

The  vertical  profiles  of  the  rms  vertical  velocity  fluctuations  in 

the  water  boundary-layer  at  fetch  x - 5 m and  for  wind  speed 

* 

• 10  m/sec. are  shown  in  Figure  37.  Data  obtained  in  three  repeated 
runs  are  Included  and,  the  repeatability  is  fairly  good.  By  comparison 
between  Figures  35  and  37,  we  obtain  that  the  ratio  of  to  v u2 

ranges  between  0,70  near  the  water  surface  and  0.50  at  larger  depths. 

The  vertical  velocity  fluctuations  again  concain  in  part  the  contribution 
from  the  orbital  motion  of  the  waves.  Figure  38  shows  the  vertical 


profiles  of  rms  and  ..he  vertical  velocity  fluctuations 


measured  at 


fetch  x ■ 5 m for  three  different  wind  speeds,  ■ 4-,  7-  and  10  m/s ec. 
Except  for  the  reduced  magnitude,  the  behavior  of  VP  is  similar  to 
that  of  Ju*  . 

The  Reynolds  shear  stress  or  the  correlation  between  the  longitudinal 
and  vertical  velocity  components  uv  is  a very  important  quantity  that 
determines  the  vertical  transfer  of  momentum  in  the  water  boundary-layer. 
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This  quantity  directly  affects  the  vertical  dispersion  of  an  oil  slick. 

Figure  39  shows  Reynolds  stress  profiles  measured  at  fetch  x ■ 5 m and 

for  three  different  wind  speeds,  U ■ 4-,  7-  and  10  m/sec.  For 

— 2 2 

U • 10  m/sec,  the  measured  maximum  value  of  -uw  is  2.6  cm  /sec 
00 

which  is  about  80  percent  of  the  shear  stress  estimated  from  the  mean 
velocity  profile  (Figure  33).  The  Reynolds  stress  has  a negative  sign 
over  most  regions  of  the  water  be undary- layer,  and,  therefore,  the 
momentum  is  transferred  downwards. 

The  turbulence  dissipation  rate  c was  determined  from  the  mean 
square  of  the  time  derivative  of  the  longitudinal  velocity  fluctuations 
by  assuming  local  isotropy  and  the  validity  of  the  Taylor’s  frozen 
turbulence  hypothesis,  i.e., 

c » 15  V (3u/3t)2/U2  . 


Figure  A0  shows  the  vertical  profiles  of  the  dissipation  rate  measured 
at  fetch  x - 5 m for  three  different  wind  speeds,  «■  4~,  7-  and 
10  m/sec. 

The  turbulence  measurements  in  the  water  boundary-layer  are  conven- 
iently tabulated  in  Tables  5 through  17  for  different  wind  speeds  and 
fetches.  Each  table  contains  the  variation  of  the  important  statistical 
quantities  with  depth  below  the  water  surface.  The  quantities  shown  are 
mean  velocity  in  the  longitudinal  direction,  root-  mean-square  velocity 
fluctuation  in  the  longitudinal  and  vertical  directions,  Reynolds  stress 
and  its  coefficient,  dissipation  rate,  dissipation  rate  weighted  by  the 
water  depth,  and  the  non-dimensional  parameter 


3/2 

(u2) 


that  has  a value  of  order  one  in  the  logarithmic  region.  It  is  noted 
that  this  parameter  becomes  larger  outside  the  boundary  layer  where  the 
signal  is  contaminated  with  instrument  noise  and  reverse  flow  effects. 
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4.7  Velocity  Spectra 


An  important  additional  statistical  quantity  that  could  be  completed 
ftom  the  random  velocity  signal  is  the  spectrum  of  the  velocity  fluctu- 
ations, or  the  distribution  of  kinetic  energy  for  different  frequency 
bands.  For  the  present  study,  such  a quantity  has  an  additional  signifi- 
cance in  that  the  contribution  of  the  wave  orbital  motion  to  the  total 
velocity  fluctuations  should  be  more  pronounced  at  the  dominant  wave 
frequency.  A typical  spectrum  S^(f)  of  longitudinal  velocity  obtained 
in  Run  No.  OIL-21  is  shown  in  Figure  41  for  fetch  x * 5 m and  wind 
speed  * 10  m/sec.  The  spectrum  shown  is  measured  at  depth  z ■ -0.8 
cm  below  the  undisturbed  surface.  The  dominant  frequency  f * 3.5  Hz, 
corresponds  exactly  to  that  in  the  wave  energy  spectrum.  It  is  seen 
from  the  figure  that  the  contribution  of  the  orbital  motion  to  the  total 
velocity  fluctuations  is  by  no  means  negligible.  The  second  harmonic 
component  at  f = 7.0  Hz  is  noticeable.  In  the  frequency  range  greater  than 
the  dominant  wave  frequency,  the  velocity  spectrum  is,  however,  wider 
than  the  wave  energy  spectrum,  an  indication  of  the  transfer  of  turbulent 
energy  from  low  frequency  (large  scale)  to  high  frequency  (amall  scale). 

The  broken  line  shown  on  the  figure  represents  the  relation 
Su(f)  * 25f  , that  is  expected  to  be  valid  for  an  inertial  subrange. 

The  -5/3  law  range  covers  about  one  decade  of  the  frequency  domain.  By 
using  Taylor's  frozen  turbulence  hypothesis  we  can  convert  the  frequency 
spectrum  Su(f)  into  the  one-dimensional  wave  number  spectrum  (k^) 

as  follows: 


U “2/3  _5/3 

E11  <kl>  " 25  %>  kl 


(4-5) 


where  k^  is  the  wave  number  component  in  the  streamwise  direction. 
Equation  (4-5)  can  then  be  used  to  estimate  the  Kolmogorov's  constant  a 
expressed  in  the  inertial  subrange  spectrum  (Hinze,  1959) 


E (k).iiae2/3k'5/3 
11  lV  55  1 


(4-6) 
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The  value  of  a is  calculated  to  be  1.79  by  inserting  l)  * 16.16  cm/sec 
2 

and  e - 109  cm  /see.  This  value  Is  reasonably  close  to  the  universal 
constant  a - 1.44  determined  from  many  other  turbulence  measurements 
(Phillips  1966).  Lin  (1974)  observed  the  universal  constant  to  be 
a • 1.66  by  fitting  experimental  data  to  his  semi-analytical  spectrum 
in  the  inertial  and  dissipation  ranges.  The  Kolmogorov  wave  number, 

* (e/v3)1^4  , is  calculated  to  be  100  cm  ^ and  the  corresponding 
Kolmogorov  length  scale  is  0.63  mm.  The  turbulence  microscale  defined 
as 


A2  - U2/(3u/3x)2 

is  A •*  0.24  cm,  thus  the  turbulence  Reynolds  number  based  on  the 
turbulence  microscale,  R^  * s/ u A/v  , is  about  150.  This  turbulence 
Reynolds  number  is  large  enough  that  the  effects  cf  the  Reynolds  number 
on  the  turbulent  motions  are  unlikely  to  be  important.  It  is  worth 
mentioning  that  in  laboratories,  researchers  have  found  that  turbulence 
follows  a "universal"  behavior  as  long  as  R^  is  larger  than  50. 

Figure  42  showB  the  longitudinal  velocity  spectra  measured  at  three 
different  depths,  z ■ -0.8,  -2.6  and  -7.6  cm,  and  at  fetch  x - 5 m,  and 
wind  speed  ■ 10  m/sec.  At  the  frequency  f - 3.5  Hz  which  corresponds 

to  the  dominant  wave  frequency,  all  these  spectra  show  peaks.  These 
peak  values  attenuate  very  rapidly  with  depth.  Thus,  at  larger  depths, 
the  contribution  from  the  orbital  motion  of  the  waves  to  the  water 
ioundary-layer  becomes  insignificant  and  the  turbulent  motion  dominates 
the  boundary  layer  flow.  Here,  we  define  turbulence  as  the  non-organized 
motion  in  contrast  to  the  organized  orbital  motion  due  to  waves.  At 
frequencies  other  than  the  dominant  wave  frequency  the  motion  is  governed 
by  the  turbulence.  The  turbulent  energy  at  frequencies  larger  than  the 
dominant  wave  frequencies  attenuates  rapidly  with  depth.  The  Kolmogorov 
frequency,  defined  as 
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is  calculated  to  be  256-,  104-  and  44  Hz  for  z ■ -0.8,  -2.6,  and  -7.6  cm  ; 

respectively.  In  other  words,  the  viscosity  cut-off  frequency  for  the  1 

I 

velocity  spectrum  decreases  with  the  increase  of  depth.  Figure  43  shows  j 

the  spectra  of  vertical  velocity  fluctuations  measured  at  x * 5 m, 

» 10  m/sec  for  three  depths,  z » -0.8,  -2.0,  and  -7,6  cm.  These 

spectra  behave,  in  general,  similarly  to  those  plotted  in  Figure  42.  ! 

i 

More  Insight  into  the  dynamics  of  the  water  boundary-layer  could  be 
gained  by  computing  the  co-  and  quad-spectra  of  the  longitudinal  and 
vertical  velocity  fluctuations.  These  are  the  real  and  imaginary  parts, 
respectively,  of  the  Fourier  transform  of  the  longitudinal  and  vertical 
velocities  cross-correlation.  An  example  of  these  spectra  is  shown  In 
Figure  44,  for  fetch  x « 5 m,  wind  speed  U ■ 10  m/sec,  and  depth 

z ■ -2.6  cm  below  the  water  surface.  It  is  seen  that  both  co-  and  quad-  ; 

spectra  peak  at  the  dominant  wave  frequency  fQ  » 3.5  Hz.  The  magnitude* 

and  phase**  of  those  spectra  are  shown  on  Figure  45.  The  magnitude 

peaks  also  at  the  dominant  wave  frequency  and  the  phase  at  this  frequency 

is  -78  deg.  This  figure  shows  that  the  shear  stress  contributes  not 

only  from  turbulent  motion,  but  also  from  the  orbital  motion  of  waves. 

These  results  indicate  the  importance  of  the  waves  in  transporting 

momentum  vertically.  ; 


Magnitude  = j (co-spectrum)  + (quad-spectrum) 


2 . , . _ . 2 j 1/2 


Phase  = tan-1  (^-spectrum 
co-spectrum 
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EFFECTS  OF  WAVES  ON  WATER  MOTION 

In  this  section,  we  present  the  experimental  results  concerning  the 
effects  of  waves  on  the  water  motion.  For  all  the  experiments,  the  wind 
speed  was  the  same,  U *■  10  m/sec,  and  waves  of  three  different  heights 
were  generated  by  the  wave  maker.  We  measured  the  wave  characteristics 
of  these  waves.  The  surface  drift  velocity  under  these  wave  conditions 
was  measured,  and  the  effect  of  wave  steepness  on  the  surface  drift 
velocity  was  determined.  The  water  boundary-layer  developed  under  the 
action  of  wind  and  waves  was  visualized  with  dye,  and  the  thickness  of 
the  water  boundary-layer  was  determined  from  the  visualization  results. 
Vertical  profiles  of  mean  velocity  and  nns  velocity  fluctuations  were 
measured  using  the  LDV.  The  velocity  spectra  were  calculated  from  the 
velocity  fluctuations  to  determine  the  wave  effect  on  the  velocity  field 
and  the  basic  turbulent  characteristics. 

5 . I Wave  Characteristics 

The  wave  maker  described  in  Section  2 was  used  to  generate  surface 
waves.  The  oscillatory  motion  of  the  wave  plunger  was  controlled  by  an 
electronic  wave  generator  operating  at  a frequency  of  about  1.6  Hz.  The 
voltages  of  the  wave  generator  were  selected  to  produce  three  different 
strokes  of  the  wave  plunger,  namely,  ■ 2.54-,  5.08-  and  7.62  cm. 
Figure  46  (a)  through  (d)  shows  the  pnotographs  of  the  water  surface  for 
Hq  - 0- , 2.54-,  5.08-  and  7.62  cm.  All  waves  propagate  from  the  left  to 
the  right.  The  photographs  were  taken  at  a fetch  of  about  x ■ 4 m.  For 
Ho  • 0 cm,  there  are  no  mechanically-generated  waves.  The  waves  are 

generated  by  wind.  For  H - 7.62  cm,  the  photograph  shows  the  spilling 

\ o 

break  up  at  the  wave  crest.  In  the  photograph,  the  breaking  region  is 
shown  as  a bright  region  across  the  tank.  The  motion  in  the  breaking 
wave  region  was  so  violent  that  we  observed  the  entrainment  of  air 
into  the  water.  Figure  47  shows  the  breaking  wave  region  viewed  from 
below  the  wave  crest.  In  this  figure,  the  direction  of  the  wave  propo- 
gation  is  from  the  left  to  the  right.  Under  the  wave  crest,  there 
appeared  bright  spots  which  indicated  air  bubbles.  These  air  bubbles 
penetrated  to  a depth  of  about  one  wave  height. 
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Figure  48  shows  Che  wave  energy  spectra  S(f)  measured  with  a 
capacitance  probe  at  fetch  x - 2 m and  wind  speed  - 10  m/sec.  The 
spectra  peak  at  the  dominant  frequency  f ■ 1.6  Hz  and  subsequent 
higher  harmonics.  The  strength  of  higher  harmonic  wave  components 
Increases  with  Hfl  . For  Hq  - 7.62  cm,  the  spectral  peaks  at  the 
dominant  wave  frequency  and  higher  harmonics  have  an  envelope  expressed 
by  S(f)  - 1.46  x 10"4  g2  f'5  for  f < f < 10  Hz. 

The  phase  speed,  Cp  , of  the  dominant  waves  was  measured  with  two 
capacitance  probes  separated  by  a distance  Ax  - 1.02  cm  apart.  These 
two  capacitance  probes  were  made  of  a thin  copper  wire  coated  with 
teflon.  We  selected  the  teflon  coated  wire  probe  instead  of  the  bare 
stainless  steel  probe  In  order  to  minimize  the  possible  cross  talk 
between  two  adjacent  probes.  To  determine  the  value  of  Cp  , we  calculated 
the  cross  correlation  of  the  surface  displacements  measured  with  the  two 
capacitance  probes.  The  time  lag  of  the  maximum  correlation.  At  , can 
therefore  be  determined,  and  it  should  correspond  to  the  time  required 
for  the  waves  traveling  over  a distance  Ax  . Hence,  the  dominant  phase 
speed  can  be  estimated  as  Cp  * Ax/ At  . The  dominant  wave  speeds  are 
108. 4-,  108.4-  and  121.1  cm/sec  for  • 2.54-,  5.08-  and  7.62  cm  respec- 

tively. Based  on  the  measurements  of  C and  f , we  can  estimate  the 

oo 

dominant  wavelength  X » Cp/fo  88  ^5. 3-  and  73.0  cm  for  the 

three  respective  values  of  Hq. 

Based  on  the  measurements  of  X , one  can  estimate  the  wave  phase 
speed  by  the  linear  wave  theory.  Comparison  between  the  direct  measure- 
ment and  the  linear  wave  prediction  could  then  indicate  the  nonlinear 
effects  on  the  wave  propagation.  Figure  49  shows  the  plot  of  Cp/  >J gX/2ti 
versus  the  wave  steepness  H/X  . Here  H is  the  wave  height  w 'ch  can 
be  estimated  from  the  rms  wave  amplitude  by  . For  comparison, 

a solid  line  representing 
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is  plotted  (Kinsman,  1965).  The  last  two  terms  in  equation  (5-1) 
represent  the  nonlinear  effect  on  the  phase  speed.  The  deviation  of  the 
measurement  points  from  the  solid  curve  could  partly  be  attributed  to 
the  wind  induced-surface  drift.  Figure  50  shows  the  variation  of  rms 
wave  amplitude  with  fetch  for  three  values  of  \\q  . We  used  the  photo- 

diode wave  gauge  to  measure  the  vertical  displacement  of  the  water 
surface  because  the  wave  height  measurement  by  the  capacitance  probe 
becomes  inaccurate  when  the  wave  breaking  occurs.  Each  data  point  was 
calculated  from  a record  of  about  60  sec.  The  overall  statistical  error 
is  estimated  to  be  about  10  percent.  This  error  is  mainly  contributed  by 
a low-frequency  modulation  (having  a period  of  about  20  sec)  of  the  wave 
field  which  is  possibly  caused  by  the  return  flow  in  the  wave  tank. 

For  Hq  ■ 2.54-  and  5.08  cm,  the  rms  wave  amplitude  Increases  with 

increasing  fetch.  For  example,  the  wa -e  energy  per  unit  density,  gn^  , 

3 2 3 2 

increases  from  980  cm  /sec  at  x * 2 m to  2509  cm  /sec  at  x ■ 5 m for 
Hq  - 2.54  cm.  In  section  4.2,  we  have  presented  the  measurement  of  rms 
wave  amplitude  for  the  case  without  mechanically-generated  waves,  but 

with  wind  only;  the  wave  energy  per  unit  density  increaoeB  from  prac- 

3 2 

tically  zero  at  x ■ 2 m to  134  cui  /sec  at  x * 5 m.  This  example  shows 
that  when  mechanically-generated  waves  are  present,  the  energy  transfer 
from  wind  to  waves  appears  to  be  more  efficient  than  when  mechanically- 


generated  waves  are  not  present.  This  observation  also  applies  to  the 

case  of  ■ 5.08  cm.  For  ■ 7.62  cm,  however,  the  wave  energy  decreases 

with  fetch.  Appsrently,  the  wa>  .»  breaking  causes  significant  energy 

dissipation  by  turbulence  whicli  overwhelms  the  energy  transfer  from  wind. 

Figure  51  shows  the  wave  energy  spectra  measured  at  two  fetches  for 

Hq  7.62  cm  and  ■ 10  m/sec.  The  wave  breaking  reduces  the  energy 

content  at  the  dominant  wave  frequency  and  increases  the  energy  at 

higher  frequencies.  For  6 Hz  < f < 40  Hz  , the  wave  energy  spectrum 

2/3  -7/3 

follows  S(f)  - 0.054  y f . The  energy  gaps  between  higher  harmonics 
are  filled  in  when  wave  breaking  occurs.  Figure  52  shows  the  wave 
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energy  spectra  measured  at  x ■ 5 tn  for  » 10  m/sec  and  for  various 

conditions  of  mechanically-generated  waves.  In  the  capillary  wave 

regime,  the  energy  increases  with  the  value  of  H . In  the  gravity  wave 

° -4  2 -5 

regime,  there  exists  an  envelope  expressed  by  S(f)  ■ 1.45  x 10  g f 
for  the  dominant  waves  and  their  higher  harmonies.  This  envelope  could, 
therefore,  represent  the  saturated  wave  spectrum  in  the  present  investi- 
gation. Figure  53  shows  the  wave  energy  spectra  measured  at  x » 5 m for 
Hq  - 7.62  cm  and  for  three  wind  speeds.  In  the  capillary  wave  regime, 
the  wave  energy  increases  with  wind  speed.  In  the  gravity  wave  regime, 
wave  energy  does  not  change  appreciably,  and  the  energy  gaps  between 
higher  harmonics  fill  in  as  wind  speed  increases. 

5. 2 Surface  Drift  Velocity 

The  surface  drift  velocity  was  measured  using  a small  paper  disc 
floating  on  the  water  surface.  The  paper  discs  were  made  the  same  way  as 
for  the  surface  drift  measurement  described  in  Section  4.3.  The  drift 
of  the  paper  discs  was  recorded  with  a 16-mm  movie  camera.  The 
standard  deviation  of  the  drift  velocity  measurements  was  typically  about 
10  to  20  percent  of  the  mean  surface  drift  velocity.  If  breaking  waves 
appear,  paper  discs  may  not  float  well  on  the  water  surface.  Frequently, 
we  did  observe  that  paper  discs  were  entrained  into  the  breaking  wave 
region  below  the  water  surface.  Consequently,  these  paper  discs  were 
discarded  when  making  the  drift  velocity  measurement.  On  the  other 
hand,  we  also  observed  that  the  standard  deviation  in  measuring  the 
surface  drift  velocity  is  large  when  breaking  waves  are  present. 

Figure  54  shows  the  mean  surface  drift  velocity  U plotted  versus 

8 

fetch  x . The  mean  surface  drift  velocity  measured  at  x ■ 2 m are 
plotted  versus  the  wave  steepness,  H/X,  in  Figure  55.  Also  plotted  is 
the  surface  drift  velocity  measured  in  the  wind  waves  at  ■ 10  m/sec. 

This  is  a limiting  case  when  the  height  of  the  mechanically-generated 
waves  approaches  zero.  The  measurement  points  show  that  at  H/X  - 0.05,  the 
surface  drift  velocity  is  lower  than  that  for  the  case  when  the  mechan- 
ically-generated waves  are  not  present.  For  H/X  > 0.05,  the  surface 
drift  velocity  increases  with  the  increase  of  the  wave  steepness. 
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The  effect  of  wave  steepness  on  the  surface  drift  has  been  Investi- 
gated by  Reiebig  et  al. , (1973).  In  their  laboratory  experiments,  the 

wind  speed  ranged  from  0-  to  7 m/sec  and  the  wave  steepness  varied  from 
0-  to  0.069.  Their  experimental  results  show  that  the  drift  mechanisms 
caused  by  wind  and  waves  ere  not  simply  additive.  At  low  wind  speeds, 
the  wave  drift  provides  an  augmentation  to  the  wind  drift,  but  at 
higher  wind  speeds,  the  waves  cause  a net  decrease  in  the  coupled  drift 
velocity.  This  wave-induced  diminishment  Increases  as  the  wind  speed 
increases.  In  the  present  investigation,  we  covered  a wider  range  of 
H/X  . Our  data  show  a similar  wave-induced  decrease  when  H/X  increases 
from  0 to  0.05.  For  H/X  > 0.05  the  surface  drift  velocity  increases 
with  wave  steepness.  In  the  following  paragraphs,  we  will  try  to  explain 
the  possible  mechanisms  which  contribute  to  the  complex  behavior  of  the 
surface  drift  due  to  the  coupled  wind  and  waves. 

Jeffreys  (1925)  hypothesized  the  separation  mechanism  for  energy 
transfer  from  wind  to  waves.  Recently,  Chang  et  al. , (1969)  conducted  a 
series  of  laboratory  measurements  of  mean  and  turbulent  velocity  in  the 
air  flow  region  above  wind-generated  waves  and  their  experimental  results 
tend  to  support  the  separation  mechanism.  Banner  and  Melville  (1976) 
conducted  a flow  visualization  experiment  to  demonstrate  that  the  air- 
flow could  separate  from  the  water  surface  in  the  lee  aide  of  mechan- 
ically-generated waves.  Based  on  kinematic  arguments,  they  hypothesized 
that  the  separating  flow  would  occur  when  waves  break.  More  recently, 
Weissman  (1978)  conducted  laboratory  experiments  in  the  airflow  over 
mechanically-generated  waves.  His  experimental  results,  both  from  smoke 
visualization  and  from  probe  measurements,  showed  that  the  separating 
flow  in  the  air  could  occur  even  though  wave  breaking  was  not  apparent. 

In  view  of  the  separation  mechanism  described  above,  we  expect  that 
in  our  experiment  conducted  for  a wind  speed  of  UM  ■ 10  m/sec  over 
mechanically-generated  waves,  the  airflow  could  separate  from  the  water 
surface  in  the  lee  side  of  the  waves.  Although  the  presence  of  the  air 
separation  could  enhance  the  energy  transfer  from  wind  to  waves,  the 
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effective  shear  stress  exerted  on  the  water  by  the  wind  is  reduced. 

This  reduction  is  not  difficult  to  imagine  if  we  realize  that  the 

separated  flow  would  introduce  a reverse  flow  region  in  which  the  shear  i 

stress  is  exerted  in  a direction  opposite  to  the  wind.  Now,  if  we 

assume  that  the  surface  drift  due  to  the  wind  should  correspond  to  the  i 

i 

effective  shear  stress  exerted  on  the  water  surface,  we  would  therefore  ! 

t 

expect  that,  for  the  case  when  mechanically-generated  waves  are  present  j 

j 

and  the  separation  exists  in  the  lee  side  of  the  waves,  the  surface  j 

drift  velocity  would  be  lower  than  for  the  case  without  mechanically-  j 

i 

generated  waves.  j 

I 

When  the  wave  steepness  increases  further,  additional  mechanisms  ! 

other  than  the  separation  mechanism  could  influence  the  surface  drift.  j 

For  example,  the  Stokes'  drift  (1847)  on  the  water  surface  is  expressed  j 

by  Phillips  (1966)  as 

i 

t 

ud  = *2  <I>2  c (5-2)  : 

where  C is  the  wave  speed  and  A is  the  wavelength.  The  Stokes'  J 

drift  would  become  important  as  the  wave  steepness  and  the  wave  speed  ! 

Increase.  Pottinger  and  Reisbig  (1973)  conducted  an  experimct.t  to  ! 

? 

determine  the  effect  of  waves  on  the  drift  of  an  oil  slick.  The  waves  ; 

were  generated  mechanically  by  a wave  maker.  The  wavelength  ranged  j 

from  25-  to  50  cm  and  the  wave  steepness  from  0.01  to  0.07.  The  oil  j 

i j 

slick  was,  in  essence,  an  oil  lens  of  diameter  4.6-  and  7.5  cm  and  thickness  i 

-4  -5  ! 

varying  from  10  to  10  of  the  wavelength.  Their  data  Indicate  that  ! 

the  oil  slick  drifted  faster  than  that  predicted  by  the  Stokes'  drift. 

They  obtained  an  empirical  relation  to  correlate  their  experimental 

results . 


U - U . + U 
s d c 


- *2(f)2  C + 0.445  (y)  C , for  0.01  < H/A  < 0.07 


(5-3) 


and  25  < A < 50  cm  , 


where  Uc  is  the  correction  term  of  the  surface  drift  velocity. 
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In  Che  following,  we  will  utilize  Che  empirical  result  expressed  in 
equation  (5-3)  to  help  us  understand  our  observation  of  the  surface 
drift  caused  by  wind  and  waves.  In  Figure  55,  we  replot  our  data  by 
subtracting  the  surface  drift  velocity  by  the  corresponding  Stokes' 
drift  velocity,  i.e.,  U - U.  , and  subsequently  by  the  correction  term 

8 a 


u 


i.e. , U 


points  of  U 


U . - U 
d c 

U . - U 
d c 


. It  is  interesting  to  point  out  that  the  data 
does  not  show  significant  variation  with  wave 


steepness  and  they  have  an  average  value  of  about  0.02  U which  corresponds 
to  wind-induced  current.  This  value  is  apparently  lower  than  that 
observed  for  the  case  when  mechanically-generated  waves  are  not  present. 
Therefore,  we  introduce  an  empirical  formula  for  the  surface  drift 
caused  by  wind  and  waves: 


U 

tr  * °-02  + 


2.H.2  , « ... 

<r  (j)  + 0.445  tj) 


C 

U 


(5-4) 


for 


0.05  < H/X  < 0.10 


10  m/sec  and  65  cm  < X < 73  cm 


A solid  line  is  drawn  in  Figure  55  to  represent  equation  (5-4)  for 
C/U,  •■0.12  , and  it  can  fairly  approximate  our  measurements. 

5. 3 Water  Boundary-Layer 

The  water  boundary-layer  was  visualized  with  dye  techniques.  The 
dye  was  introduced  into  the  water  by  a vertical  string  0.3  cm  in  diameter. 
The  string  was  located  along  the  centerline  of  the  tank.  The  string  was 
made  of  cotton  and  formed  a closed  loop  with  three  point 
supports.  Part  of  the  string  was  stretched  vertically  in  the  water  and 
part  of  it  was  exposed  in  the  air.  The  part  of  string  exposed  in  the 
air  was  dyed  with  a fluorescent  dye,  and  it  was  then  lowered  below  the 
water  by  sliding  the  string  loop.  The  dye,  washed  away  from  the  string, 
was  initially  vertical,  and  this  vertical  dye  line  was  deformed  and 
dispersed  by  the  motion  in  the  water.  Figure  56  is  a photograph  showing 
the  dye  dispersion  in  the  water  for  Hq  - 7.62  cm  and  *>  10  m/sec. 

Two  flow  regions  are  distinctly  visible.  In  the  top  region  about  15  cm 
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below  the  water  surface,  Che  dye  disperses  it  the  wind  direction,  and  in 
Che  bottom  region  below  15  cm  from  the  water  surface  Che  dye  disperses 
in  the  reverse  direction.  The  top  region  is  apparently  the  boundary 
layer  of  flow  set  up  by  wind  and  waves,  and  the  bottom  region  corresponds 
to  the  return  flow  due  to  the  finite  length  of  the  tank.  For  the  case 
when  waves  are  generated  by  wind  only,  we  used  a false  bottom  to  channel 
the  return  flow  in-between  the  tank  bottom  and  the  false  bottom,  and 
therefore  we  were  able  to  control  and  minimize  the  reverse  flow  in  the 
region  between  the  water  surface  and  the  false  bottom.  We  tried  to  use 
the  same  setup  of  the  false  bottom  for  the  case  when  both  wind  and 
mechanically-generated  waves  were  present,  but  the  reverse  flow  could 
not  be  controlled  as  well  a-  when  only  wind-generated  waves  were  present. 
For  lack  of  time  to  try  to  find  a better  way  to  control  the  reverse 
flow,  we  decided  not  to  use  the  false  bottom  for  the  case  when  both  wind 
and  mechanically-generated  waves  were  present. 

The  presence  of  the  reverse  flow  could  affect  the  water  boundary- 
layer.  We  are  not  in  a position  to  estimate  these  effects  due  to  the 
reverse  flow  quantitatively,  but,  on  the  other  hand,  we  expect  that  in 
the  region  near  the  water  surface  and  away  from  the  reverse  flow,  the 
flow  characteristics  should  be  most  dominated  by  the  local  flow  conditions 
such  as  the  break  up  of  waves  and  the  presence  of  strong  shear.  Figure 
57  presents  the  thickness  of  the  water  boundary-layer  measured  from  the 
movie  record  of  the  dye  dispersion  in  the  water  boundary-layer.  We 
define  the  thickness  of  the  water  boundary-layer  as  the  depth  measured 
between  the  undisturbed  water  surface  and  the  vertical  position  where 
the  local  mean  velocity  was  approximately  zero.  The  data  do  not  show 
any  consistent  varistion  with  fetch  and  wave  height. 

5.4  Mean  Velocity 

Figures  58  and  59  present  the  vertical  profiles  of  longitudinal  and 
vertical  mean  velocities  for  Hq  ■ 7.62-  and  2.54  cm.  In  general,  we 
experienced  some  difficulties  in  measuring  mean  velocity  components 
because  of  the  unsteady  nature  of  the  return  flow  set  up  in  the  tank. 
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This  unsteady  flow  apparently  does  not  affect  the  measurement  of  velocity 
fluctuations  (Section  5.6),  but  it  docs  affect  the  mean  velocity  component. 
This  unsteady  phenomenon  has  a typical  period  of  about  20  sec,  but  our 
measurement  period  was  about  40  sec.  Thus,  we  do  experience  some 
statistical  uncertainty  in  measuring  the  mean  velocity.  Figure  58  shows 
that  the  mean  velocity  profile  has  a drift  layer  in  the  region  about 
15  cm  below  the  undisturbed  water  surface,  and  a reverse  flrv  between 
the  depths  of  15-  and  25  cm.  The  drift  layer  does  not  appear  as  strong 
as  when  the  mechanically-generated  waves  are  absent  (Figure  32).  This 
result  shows  that  in  the  presence  of  mechanically-generated  waves,  wind 
is  no  longer  as  effective  in  transferring  energy  into  current  as  when 
the  mechanically-generated  waves  are  absent.  On  the  other  hand,  we  have 
observed  that  when  mechanically-generated  waves  are  present  the  energy 
transfer  from  wind  to  waves  appears  to  be  more  efficient  than  when 
mechanically-generated  waves  are  absent  (Section  5.1).  A systematic 
measurement  of  both  energetics  in  the  air,  water  and  wave  fields  is 
needed  to  clarify  how  the  energy  is  transferred  from  the  air  to  the 
drift  current  and  to  waves. 

Figure  59  shows  that  in  the  top  drift  layer  the  vertical  mean 
velocity  is  small,  while  in  Che  reverse  flow  region  the  vertical  velocity 
becomes  significant  when  compared  with  the  longitudinal  velocity. 

5.5  RMS  Velocity  Fluctuations 

Figure  60  shows  the  vertical  profile  of  rms  longitudinal  velocity 

fluctuations,  -J  measured  for  H » 7.62-  and  2.54  cm.  For 

o 

Hq  * 7.62  cm,  the  data  points  are  denoted  by  circles.  At  several  depths, 
repeated  measurements  were  made  to  examine  the  repeatability  in  the 
measurement  of  rms  velocity  fluctuations.  For  Hq  ■ 2.54  cm,  repeated 
measurements  were  also  made  in  two  series  of  experiments.  This  good 
repeatability  shows  chat  the  measurement  of  rms  velocity  fluctuations 
does  not  appear  to  be  significantly  affected  by  the  return  flow. 

Especially  in  the  drift  layer,  the  reverse  flow  should  not  appreciably 
influence  the  velocity  field.  Comparison  between  Figures  58  and  60 
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shows  that  the  rms  velocity  fluctuations  are  much  larger  than  the  mean 
velocity  components.  This  situation  is  in  contrast  with  the  case  when 
the  mechanically-generated  waves  were  not  present;  the  mean  velocity 
components  are  larger  than  the  corresponding  rms  velocity  fluctuations 
(Figures  32  and  35).  Therefore,  we  have  established  the  experimental 
condition  under  which  waves  dominate  wind  in  the  generation  of  water 
motion . 

The  vertical  profiles  of  rms  vertical  velocity  fluctuations  are 
plotted  in  Figure  61.  The  rms  vertical  velocity  fluctuations  behave  in 
a similar  manner  as  the  rms  longitudinal  velocity  fluctuations,  and 
their  magnitudes  are  about  the  same.  Figure  62  is  the  replot  of  Figures 
60  and  61  on  a seni-logarithmic  paper.  The  measurements  can  be 
approximated  by 

J *=  33  exp  (-|z|/'10.9),  for  » 7.62  cm  and  3.5  < -z  < 17  cm 

(5-5) 

and 

V ■ 15  exp  (—  j z | / 10. 9> , for  Hq  ■ 2.54  cm,  2.5  < -z  < 20  cm. 

(5-6) 

This  exponential  attenuation  of  ; 3 velocity  fluctuations  demonstrates 
clearly  that  the  velocity  field  is  dominated  by  wave  motion.  The  constant 
10.9  expressed  in  equations  (5-5)  and  (5-6)  represents  the  characteristic 
depth  of  the  velocity  attenuation  and  it  is  about  1/6  of  the  dominant 
wave  length.  At  depths  larger  than  about  20  cm  or  about  1/4  to  1/3  of 
the  dominant  wavelength,  the  rms  values  attenuate  at  a rate  slower  than 
the  exponential  rate. 
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5.6  Mean  Square  Velocity  Time  Derivatives 

Figure  63  shows  vertical  profiles  of  the  mean  uqitare  velocity 

time  derivatives.  In  a stationary  flow  in  which  the  mean  velocity  is 

much  larger  than  the  fluctuating  velocity,  the  mean  square  velocity  time 

derivatives  (3u/&t)^.  can  be  related  to  (3u/3x)2  according  to  Taylor’s 

frozen  turbulence  hypothesis,  i.e.,  (3u/3t)2  = (3u/3x)^  . Hence,  by 

further  assuming  local  isotropy,  we  can  estimate  the  turbulent  dissipation 

rate.  For  the  present  flow  condition,  the  rms  velocity  fluctuation  is 

larger  than  the  mean  velocity,  therefore  the  application  of  the 

Taylor's  hypothesis  is  questionable.  It  is  our  belief,  however,  that 

the  measurement  of  <8u/3t)^  or  (3w/3c)2  could  still  indicate  small 

scale  turbulence  motion.  Especially  for  the  present  situation,  the 

velocity  fluctuation  is  dominated  by  the  wave  motion,  thus  the  rms 

values  of  velocity  fluctuation  cannot  be  used  to  indicate  turbulent 

motion.  Although  (3u/3t)2  or  (3w/3t)2  cannot  directly  relate  to  the 

dissipation  rate,  its  measurement  can  give  the  local  acceleration. 

In  case  H *2.54  cm,  the  roaximum\/  (3w/3t)2  is  about  141  ci.,/sec2  which 
o 

is  about  0.14  times  the  gravitational  acceleration.  But,  for  H - 7.62  cm, 
che  maximum^  (3w/3t)^  becomes  as  high  as  1680  cm/sec  which  is  about 
twice  the  gravitational  acceleration.  In  general,  we  observed  that  when 
wave  breaking  occurs,  extremely  high  dissipation  rate  or  large  local 
fluid  acceleration  appear. 


5.  7 Velocity  Spectra 

Figure  64  shows  the  spectrum  of  longitudinal  velocity  fluctuations 
measured  at  z * - 3.5  cm.  Two  distinct  spectral  regions  exist.  Near 
the  dominant  wave  frequency  f ■ 1.6  Hz,  the  spectrum  peaks  and  is 
dominated  by  the  wave  motion.  At  frequencies  higher  than  the  dominant 
wave  frequency,  the  spectrum  is  distributed  over  a wide  frequency  range 
therefore,  it  represents  a turbulence  motion.  In  particular,  the  -5/3 
power  law  can  be  fitted  as 


Su(f)  - 48  f"5/3,  for  10  < f < 100  Hz  . 


(5-7) 


In  Figure  65,  the  spectrum  of  vertical  velocity  fluctuations  is 
shown.  It  has  a dominant  peak  and  a second  harmonic  component  which 
does  not  appear  in  the  spectrum  of  longitudinal  velocity  fluctuations. 

For  2 < f < 10,  the  spectral  density  of  the  vertical  velocity  fluctu- 
ations is  higher  than  that  of  the  horizontal  velocity  fluctuations.  At 
high  frequencies,  the  spectrum  in  Figure  65  is  fitted  by 

S (f)  - 66  f~5/3,  for  4 < f < 100  Hz  . (5-8) 

w 

Comparison  between  Figures  64  and  65  indicates  a possible  generation 
mechanism  of  turbulence  in  breaking  waves.  As  will  be  shown  in  Section 
6,  the  wave  breaking  introduces  very  strong  vertical  motion  as  a plume 
injected  downwards  from  the  point  where  the  wave  breaking  occurs.  The 
spectrum  of  this  vertical  motion  corresponds  to  the  energy  excess  in  the 
spectrum  of  vertical  velocity  fluctuations  for  2 < f < 10.  The  spectrum 
shown  in  Figure  65  is  identical  to  that  shown  in  Figure  64  for  the 
frequency  range  1 < f < 2 Hz  . The  identical  spectra  for  S^tf)  and 
Sw(f)  shows  unmistakably  that  in  the  frequency  range  1 < f 5 2 Hz 
the  motion  is  caused  by  waves. 

In  the  frequency  range  2 _<  f < 4 Hz,  Su(f)  *8  l°wer  than  S^ff), 

and  the  slope  of  S (f)  is  smaller  than  that  of  S (f).  Based  on  these 
u w 

observations,  we  speculate  that  strong  energy  transfers  from  the  vertical 
velocity  component  to  the  longitudinal  velocity  component.  This  spec- 
ulation can  readily  be  proven  or  disproven  if  the  co-  and  quad-spectrum 
of  longitudinal  and  vertical  velocity  fluctuations  car.  simultaneously  be 
measured. 

From  equations  (5-7)  and  (5-8),  the  ratio  Sw(f)/S  (f)  is  estimated 
to  be  about  1.38,  This  value  is  very  close  to  the  value  of  4/3  when  the 
turbulence  is  locally  isotropic.  Although  we  question  the  applicability 
of  Taylor's  hypothesis  to  the  turbulence  in  breaking  waves,  it  is  still 
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very  convenient  to  use  this  hypothesis  so  that  we  can  estimate  the 
dissipation  rate.  Two  methods  are  suggested: 

(i)  By  fitting  the  measured  S^ff)  to  the  Kolmogorov  spectrum 

(Equation  4-6).  Let  U ■ 7.4  cm/ sec,  a ■ 1.66,  then 

- 2 3 

e ■ 738  cm  /sec  . 

(ii)  By  using  the  formula 

— i (— )^  _ 2 3 

c * 15v  —2  v3t  , then  e = 4629  cm  /sec  . 

In  view  of  the  uncertainties  in  the  applicability  of  Taylor's  hypothesis, 

the  values  of  e estimated  by  these  two  methods  are  reasonably  accepted. 

From  the  above  estimates  of  e , one  can  calculate  the  Kolmogorov  length 

scale  to  be  0.06  or  0.04  mm,  respectively. 

The  Taylor  microscale  is  estit.  ted  to  be  about  0.1  cm  and  the 

corresponding  turbulent  Reynolds  number,  R^,  is  about  200.  This  value 

of  is  so  high  that  the  viscosity  dissipation  is  negligible  when 

compared  with  turbulent  dissipation.  Hence,  the  effect  of  Reynolds 

number  is  not  important  in  the  present  investigation. 

Figures  66  and  67  present  the  spectra  of  longitudinal  and  vertical 

velocity  fluctuations  measured  at  z ■ -10  cm.  Both  spectra  peak  at  the 

-5/3 

dominant  wave  frequency  and  have  turbulence  spectra  ® f 
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6. 


viIL  DISPERSION 


Based  on  che  measurements  of  wave  characteristics,  surface  drift 

velocity,  and  velocity  in  the  water  boundary-layer  for  the  case  when 

both  wind  and  mechanically-generated  waves  are  present  we  selected  the 

experimental  conditions  Uw  " 10  m/sec  and  Hq  * 7.62  cm  for  the  oil 

dispersion  experiments.  Diesel  oil  No.  2 was  usi  for  the  experiments. 

The  oil  discharge  was  about  50  cc/sec.  The  specific  gravity  was  0.84, 

2 

che  kinematic  viscosity  was  0.19  cm  /sec.  The  thickness  of  the  oil 
slick  at  the  discharge  point  was  not  measured,  but  it  was  estimated  as 
follows.  By  assuming  that  the  oil  slick  drifts  at  the  surface  drift 
velocity  (Figure  54)  and  that  the  slick  is  uniform  across  the  tank,  we 
estimate  the  slick  thickness  to  be  about  0.1  mm.  In  reality,  the  oil 
slick  thickness  was  not  distributed  uniformly  because  of  the  presence  of 
the  boundary  layer  of  the  oil  slick  by  the  side  walls.  We  also  observed 
that  the  oil  slick  drifts  in  the  form  of  streaks  which  resemble  the  wind 
rows  observed  in  wind-generated  waves. 

6 . 1 Flow  Visualization 

Flow  visualization  experiments  were  conducted  to  investigate  the 
evolution  of  the  oil  dispersion  as  the  oil  slick  drifted  downwind.  The 
visualizaton  results  were  recorded  using  a 35-mm  camera  and  a 16-mm 
movie  camera.  Initially, ' the  oil  alick  was  in  sheet  form.  As  the  oil 
slick  drifted  downwind  under  the  action  of  wind  and  breaking  waves,  the 
oil  slicks  were  broker,  into  patches  and  lenses.  Because  of  the  combined 
effects  of  gravity  and  surface  tension,  the  shape  of  the  oil  lenses  was 
approximately  elliptical  with  ’arger  horizontal  dimension  than  vertical 
dimension.  Furtiier  downstream,  oil  droplets  of  smaller  sizes  were 
formed  either  by  the  break  up  of  the  oil  lenses  or  by  the  break  up  of 
the  thin  oil  patches  (sheet-like) . These  droplets  were  spherical  and 
they  were  present  in  bursts  at  a frequency  corresponding  to  the  dominant 
frequency  of  the  breaking  waves.  The  burst  of  oil  droplets  appeared  as 
a plume  plunging  downwards  below  the  water  surface.  The  horizontal 
extent  of  the  oil  droplet  burst  decreased  with  depLh.  Host  of  the  oil 
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droplets  resurfaced  very  rapidly,  and  oil  droplets  of  small  sizes  stayed 
in  the  water  boundary-layer  and  convected  downwind  by  the  drift  currents. 

From  time  to  time,  we  observed  droplets  having  sizes  of  about  0.5- 
to  1 cm  drifting  very  freely  in  the  water  as  if  they  were  nearly 
neutra’ly  buoyant.  We  speculated  that  these  droplets  were,  in  essence, 
a water  droplet  coated  with  a thin  oil  film.  We  further  observed  that 
the  oil  film  on  the  water  droplet  was  not  uniform  in  thickness.  Near 
the  top  of  the  droplet,  the  oil  film  appeared  thicker  than  that  near  the 
bottom  due  to  the  buoyancy  of  the  oil  in  water.  We  called  these  droplets 
"oil  bubbles"  because  of  their  resemblance  to  soap  bubbles  in  the  air. 

It  must  be  reminded  that  the  soap  bubble  has  a thicker  film  near  the 
bottom  of  the  bubble.  Sometimes  we  observed  that  these  "oil  bubbles" 
were  attached  to  one  or  several  air  bubbles.  These  air  bubbles  were 
entrained  into  water  by  the  breaking  waves.  The  conglomerate  of  the  oil 
and  air  bubble  resurfaces  very  rapidly  because  of  t'ne  additional  buoyancy 
of  the  air  bubble. 

6 . 2 Oil  Concentration 

We  used  the  photographic  technique  to  measure  the  vertical  profile 
of  oil  concentration  in  the  water  boundary-layer.  The  sample  volume  for 
oil  concentration  measurements  was  selected  at  fetch  x “ 4 m.  The 
sample  volume  was  defined  between  the  side  (glass)  wall  and  a black 
colored  vertical  plate.  The  plate  was  about  30  cm  high,  2 cm  long  and 
0.2  cm  thick  and  was  positioned  1.65  cm  away  from  and  parallel  to  the 
side  wall.  The  total  sample  volume  was  a rectangular  column,  30  cm 
high,  2 cm  long  and  1.65  cm  thick.  Flood  lamps  projected  light  from 
above  for  illumination.  The  oil  droplets  existing  in  this  sample  volume 
were  recorded  by  a polaroid  camera  standing  about  30  cm  away  from  the 
side  of  the  tank.  To  determine  the  volume  occupied  by  the  oil  droplets, 
their  size  was  meusured  and  their  number  was  counted  under  a microscope. 
The  volume  measurements  were  made  at  vertical  increments  of  1.25  cm.  At 
each  increment,  the  oil  concentration  wa3  determined  by  dividing  the  oil 
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volume  to  the  total  volume  of  the  increment,  4.1  cm  . The  oil  droplets 
observed  in  this  figure  are  approximately  spherical  in  shape,  thus  by 
measuring  the  diameters  (size)  of  the  oil  droplets,  one  can  calculate 
the  volume  accordingly.  Identifying  the  oil  droplets  in  this  figure  is 
by  no  means  a trivial  task.  Under  these  light  conditions,  an  oil  droplet 
is  accompanied  by  two  bright  spots,  one  on  the  top  and  one  on  the 
bottom  of  the  droplet.  The  top  bright  spot  results  from  the  light 
reflection  from  the  top  surface  of  the  oil  droplet,  and  the  bottom 
bright  spot  results  from  the  light  passing  through  and  refracting  out  of 
the  oil  droplet. 

The  vertical  profile  of  instantaneous  concentration  is  plotted  in 
Figure  68.  This  vertical  profile  shows  that  the  maximum  volume  concen- 
tration Is  about  0.009  and  that  the  oil  can  penetrate  as  deep  as  about 

20  cm  (2.6  H ).  To  obtain  statistics  of  oil  concentration,  such  as  mean 
o 

and  rms,  we  used  seven  vertical  profiles  of  instantaneous  oil  concen- 
tration. These  profiles  were  measured  from  photographs  which  were 
sequentially  taken  at  time  increments  of  about  10  sec.  Hence,  each 
photograph  could  be  treated  as  an  independent  observation.  Figure  69 
shows  the  vertical  profiles  of  mean  and  rma  oil  volume  concentration,  C 
and  •J  C2  . The  rms  oil  concentration  was  measured  by  calculating  the 
standard  deviation  of  the  seven  instantaneous  concentration  profiles, 
thus  it  includes  both  wave-induced  and  turbulent  fluctuations.  The  rms 
oil  concentration  is  comparable  to  the  mean  value. 

The  vertical  profile  of  the  mean  oil  concentration  is  replotted  on 
a semi-logarithmic  graph  in  Figure  70.  The  profile  can  be  approximated 
by 


C = 0.015  exp  (-|z|/7.9),  for  5 cm  < -z  < 20  cm.  (6-1) 

The  constant  in  the  exponent  of  equation  (6-1)  is  a characteristic 
depth  7.9  cm  for  the  oil  concentration.  This  characteristic  depth  is 
comparable  to  that  observed  from  the  vertical  profile  of  rms  velocity 
fluctuations  (Section  5.5). 
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From  the  seven  photographs  aforementioned,  we  obtained  the  size 
distribution  of  oil  droplets  for  - 10  ra/sec  and  x ■ 4 m.  For  each 
photograph,  we  counted  the  total  number  of  oil  droplets.  The  oil 
droplets  were  thus  grouped  according  to  their  eizes.  The  smallest  size 
of  oil  droplets  which  could  be  resolved  photographically  was  0.5  mm,  and 
the  size  increment  was  0.25  nnn.  For  each  size  group,  we  obtained  the 
number  of  oil  droplets.  By  dividing  this  number  by  the  total  number  of 
oil  droplets,  we  can  then  obtain  the  percentage  frequency  for  that  size 
of  oil  droplets;  therefore,  the  size  distribution  of  oil  droplets  for 
this  photograph  is  determined.  From  the  seven  photographs,  we  can 
further  obtain  the  average  percentage  frequency  for  a specific  size. 
Figure  71  is  the  size  distribution  of  oil  droplets.  The  average  oil 
droplet  size  is  about  1.1  mm  and  the  most  frequently  observed  size  is 
0.75  mm  for  the  present  investigation. 
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7.  CONCLUSIONS 


7.  CONCLUSIONS 

7.1  Turbulence  and  Wave-Induced  Motion 

Based  on  the  experimental  data  obtained  in  the  present  research 
study,  we  conclude  that  the  turbulence  structure  in  surface  waters 
depends  strongly  on  its  generating  mechanisms.  For  the  present  interest, 
two  mechanisms,  namely  wind  and  waves,  were  considered.  The  important 
findings  are  summarized  in  the  following: 

(i)  When  wind  dominates  the  turbulence  generation,  the  water 

boundary-layer  behaves  very  similar  to  the  turbulent  boundary- 
layer  over  a rough  surface.  The  mean  velocity  has  a logarithmic 


region.  The  water  boundary-layer  thickness 
fetch  and  is  approximated  by 


grows  with 


~ - 0.025  (^f)475  , for  0.2  < & < 3 . 

U U U 

CD  CO  00 

The  profile  of  turbulence  intensity  also  resembles  that  which 

na9  been  observed  in  the  turbulent  boundary-layer  over  a rough 

surface.  The  maximum  turbulence  Intensity  is  about  0.7  percent 

of  wind  speed  or  20  percent  of  the  surface  drift  velocity. 

The  contribution  of  wave  motion  to  the  turbulence  intensity 

measurement  is  not  negligible.  The  velocity  spectrum  shows  a 

-5/3 

peak  at  the  dominant  wave  period  and  a modest  range  of  f law. 
(ii)  When  waves  dominate  the  turbulence  generation,  the  water  flow 
does  not  resemble  the  boundary-layer  flow.  The  water  motion 
ia  mainly  governed  by  waves.  The  wave  induced  motion,  such  as 
the  rms  velocity  fluctuation  attenuates  exponentially  with 
depth.  The  turbulent  motion  is  predominately  introduced  by 
the  wave  breaking.  The  wave  breaking  introduces  a strong 
vertical  burst  motion  like  a plume  plunging  downwards  below 
the  water  surface.  The  penetration  of  tins  burst  motion  is 
about  one  wave  height.  The  spectrum  show  a dominant  peak 
representing  the  wave-induced  motion  and  an  extensive 


Kolmogorov  spectrum  f 


for  almost  two  decades  of  frequencies. 
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The  mean-square  velocity  time  derivatives  indicate  significant 
turbulence  generation  near  the  water  surface  when  breaking 


waves  occur. 

7.2  Wave  Growth  Characteristics 

The  wave  growth  characteristics  have  been  investigated  by  measuring 
the  variation  of  the  rms  wave  height  with  fetch  for  various  wind  and 
wave  conditions.  The  important  findings- are  summarized: 

(i)  For  wind  waves,  the  present  results  are  consistent  with 
previous  measurements. 

(ii)  For  mechanically-generated  waves  with  wind,  the  waves  grow  at 
a faster  rate  than  without  mechanically-generated  waves  but 
with  wind.  In  other  words,  in  the  presence  of  mechanically- 
generated  waves,  the  energy  is  more  effectively  transferred 
from  wind  to  waves.  When  wave  breaking  occurs,  we  observe  a 
decrease  of  wave  energy  with  fetch.  Apparently,  the  dissipation 
mechanism  due  to  wave  breaking  overwhelms  the  energy  transferred 
from  wind  to  waves. 


7. 3 Surface  Drift  Velocity 

Paper  discs  were  used  as  floats  to  measure  the  surface  drift 
velocity.  The  important  findings  are  summarized: 

(i)  For  wind  waves  case,  the  present  result  is  consistent  with 

previous  results,  namely,  the  surface  drift  velocity  is  about 
3.2  percent  of  the  wind  speed. 

(il)  For  mechanically-generated  waves  with  wind,  the  surface  drift 
velocity  not  only  varies  with  wind  speed,  but  also  with  wave 
steepness.  We  obtain  an  empirical  formula  under  the  present 


experimental  conditions,  i.e 
U 


~ - 0.02  + 


tt2  (y) 2 + 0.445  (|)|  Jr 


for  0.05  < H/X  < 0.10,  V “ 10  m/sec  and  65  cm  < X < 73  cm 
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7.4  011  Dispersion 

Oil  dispersion  was  investigated  under  conditions  when  breaking 
waves  occur.  As  the  waves  break,  strong  vertical  motion  Is  Introduced 
and  the  oil  slick  is  broken  into  oil  droplets.  We  measured  the  oil 
concentration  by  using  the  photographic  technique.  For  the  particular 
parameters  used  in  the  present  experiment , the  mean  oil  concentration 
decreases  exponentially  as 

C » 0.015  exp  (— 1 2 1 / 7 . 9) » for  5 cm  < - 2 < 20  cm  , 

U * 10  m/sec  and  H ■ 7.62  cm  . 

OD  O 

The  nns  oil  concentration  has  values  comparable  to  the  mean  oil  concen- 
tration. The  size  distribution  of  oil  droplets  was  determined  photo- 
graphically. The  average  si2e  is  about  1.1  mm. 


47 


1 

•a 


1 

i 


4 

1 


3 


| 

1 


A 


8.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


Before  we  Initiated  this  research  study,  very  little  was  known  con- 
cerning the  turbulent  motion  in  surface  waters  and  the  effect  of  tur- 
bulence on  oil  dispersion.  The  present  study. provides  some  important 
Information  necessary  for  understanding  the  turbulent  dispersion  of  oil 
slicks,  and  it  also  calls  for  additional  research  effort. 

8. 1 Surface  Drift  V'elocity 

In  Section  5.3,  we  have  shown  the  complicated  coupling  effects  of 
wind  and  waves  on  the  surface  drift  velocity.  Also,  we  have  shown  that 
it  is  possible,  empirically,  to  separate  the  effects  of  wind  and  waves 
as  expressed  in  equation  (5-4).  But,  the  expression  is  only  for  a 
modest  range  of  wave  steepness  under  given  winu  speed  and  wavelength. 

It  would,  therefore,  be  interesting  to  obtain,  if  possible,  a similar 
empirical  formula  for  a variety  of  wind  and  wave  conditions.  This 
formula  could  be  used  in  numerical  modeling  for  predicting  the 
horizontal  spread  of  on  oil  slick  due  to  the  surface  dr*ft. 

8 . 2 Turbulent  Mot-ion  Versus  Wave-Induced  Motion 

Measurements  of  both  longitudinal  and  vertical  velocity  components 
in  wind  waves  and  in  mechanically-generated  waves  with  wind  show  two 
basic  motions,  namely,  turbulent  motion  and  wave-induced  motion.  The 
interaction  between  the  turbulent  motion  and  the  wave-induced  motion  is 
not  well  understood.  On  the  other  hand,  the  updated  turbulent  theory 
cannot  include  the  wave-induced  motion.  The  following  approaches  are 
suggested : 

(i)  Separate  the  turbulent  motion  from  the  wave-induced  motion  by 
the  phase-averaging  technique,  and  then  investigate  the 
characteristics  separately  for  either  type  of  motion. 

(ii)  extend  the  present  turbulence  theory  by  including  the  wave- 
induced  motion.  For  the  study  of  oil  dispersion,  we  do  not 
recommend  a complicated  turbulent  theory;  instead,  a theory 
equivalent  to  the  eddy  viscosity  approximation  could  be 
considered. 
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8.3  Taylor's  Frozen  Turbulence  Hypothesis 

Most  of  turbulence  measurements  was  done  In  time  series.  Taylor's 
frozen  turbulence  hypothesis  is  frequently  used  to  convert  statistics 
measured  In  the  time  domain  Into  the  spatial  domain.  For  the  case  when 
the  mean  velocity  is  ouch  larger  then  the  fluctuating  component, 

Taylor's  hypothesis  is  valid.  But,  when  the  fluctuating  velocity 
component  is  larger  than  the  mean  velocity  component,  the  validity  of 
Taylor's  hypothesis  is  subject  to  doubt.  Theoretical  and  experimental 
studies  are  needed  to  investigate  the  validity  of  Taylor's  hypothesis. 

8.4  Oil  Concentration  Flux 

Visualization  techniques  have  been  used  to  measure  oil  concen- 
tration. The  technique  can  further  be  uaed  to  determine  the  oil  concen- 
tration fluxes  simply  by  multiplying  the  volume  of  the  oil  droplet  by 
its  corresponding  velocity  components.  The  oil  concentration  fluxes 
could  yield  the  important  information  such  as  eddy  diffuaivity  of  oil 
dispersion. 
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APPENDIX  A.  ANNOTATIONS 


Aitken  (1884)  studied  the  effect  of  an  oil  film  on  the  waves  of  the 
sea.  He  was  motivated  by  an  earlier  experiment)  conducted  in  the 
Aberdeen  Harbour,  which  showed  the  smoothing  effects  of  oil  in  prevent- 
ing the  formation  of  dangerous  crests. 

Jeffreys  (1925,  1926)  studied  the  mechanisms  of  air-water  interac- 
tion. He  attempted  to  explain  the  physical  process  involved  in  the 
generation  of  waves  by  wind  by  using  the  separation  mechanism  hypothe- 
sis. 

Keulegan  (1951)  measured  the  surface  drift  current  in  a wind-wave 

tank. 

Longuet-Hipglna  (1953),  in  his  experiments,  indicated  that  the 
mass-transport  velocity  can  be  very  different  from  that  determined  by 
Stokes  on  the  assumption  of  a perfect  fluid.  Taking  into  account  the 
viscosity  and  boundary  conditions,  he  derived  a formula  for  calculating 
the  mass  transport  velocity  both  for  standing  and  progressive  sinusoidal 
waves . 

Miles  (1957,  1959a,  1959b,  1962,  1967)  theoretically  examined  the 
momentum  transfer  and  wave  growth  for  a sheared  airflow.  He  posed  the 
problem  as  a hydrodynamic  instability  with  small  perturbations  and  he 
obtained  a prediction  of  the  pressure  field  influencing  the  water 
surface. 

Grant,  et  al.  (1963)  studied  the  important  mechanisms  for  the 
generation  of  turbulence  in  surface  water.  He  Indicated  that  the  tur- 
bulent dissipation  rate  in  the  region  near  the  surface  is  larger  than 
that  some  distance  below  the  surface. 

Glover  (1964)  presented  a solution  for  the  diffusion  equation  in 
open  channels  with  the  initial  condition  of  an  instantaneous  point 
source  at  the  origin. 

Baines  and  Knapp  (1965)  measured  the  mean  drift  currents  induced  by 
wind  in  a shallow  channel. 
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Bye  (1965)  studied  wind-driven  circulation  in  unatratified  lakes. 

Colonell  (1966)  showed  experimentally  that  Miles's  theory  for  the 
momentum  tranpfer  and  wave  growth  of  a sheared  airflow  under-predlcts 
the  growth  rates  of  the  waves. 

Hldy  and  Plate  (1966)  observed  the  effects  of  wind  on  water  standing 
in  a laboratory  channel. 

Phillips  (1966)  reviewed  state-of-the-art  knowledge  on  wind  waves 
and  the  dynamics  of  the  upper  ocean. 

Sayre  and  Chang  (1966)  reported  a laboratory  experiment  on  the 
dispersion  process  for  dissolved,  suspended,  and  floating  dispersants  in 
an  open  channel.  They  adopted  Aria's  moment  distribution  method  to 
solve  the  diffusion  equation  numerically. 

Bole  and  Hsu  (1967)  conducted  a laboratory  experiment  to  determine 
the  response  of  waves  to  wind  excitation.  They  observed  higher  growth 
rates  for  the  waves  than  those  predicted  by  Miles's  theory. 

Fischer  (1967,  1969)  studied  the  mechanics  of  dispersion  in  natural 
streams  and  the  effect  of  bends  cn  dispersion.  He  estimated  the  depth- 
averaged  lateral  diffusion  coefficient  in  curved  channels,  and  he  derived 
an  equation  for  predicting  the  longitudinal  dispersion  coefficient  in 
wide  open  channels  by  adopting  Taylor's  convective  dispersion  theory. 

Davies,  et  al.  (1968)  studied  the  ripple  formation  caused  by  low- 
velocity  air  flowing  over  water. 

Karaki  and  Hsu  (1968)  conducted  an  experiment  to  investigate  the 
structure  of  a turbulent-wind  boundary-layer  over  a water  surface, 

Shontlng  (1968)  measured  the  vertical  velocity  fluctuations  in  wind 
waves.  The  rms  value  of  the  vertical  velocity  fluctuations  decreases 
exponentially  with  depth  for  the  water  region  having  depth  less  than  6 
m.  The  autospectrum  of  the  vertical  velocity  fluctuations  had  a dominant 
peak  at  a frequency  corresponding  to  the  spectral  peak  of  the  wave  energy 
spectrum. 
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Wu  (1968,  1969,  1975)  in  a aeries  of  reports  studied  wind-wave 
interaction  in  a wind-wave  tank.  He  measured  the  profile  of  the  drift 
current  induced  by  wind  and  estimated  the  effect  of  the  wind  on  disper- 
sion. He  pointed  out  that  waves  do  have  a strong  effect  on  the  diffusion 
process. 

Canevari  (1969)  studied  the  mechanisms  of  and  the  basis  for  dis- 
persing oil  slicks  with  chemical  dispersants. 

Chang  (1969)  investigated  the  detailed  structure  of  airflow  above 
and  between  crests  of  progressing  wind  waves. 

Eloubaldy  (1969)  studied  the  diffusion  process  in  open  channels. 

He  noted  that  waves  have  a significant  effect  on  the  diffusion  process. 

Fay  (1969,  1971)  presented,  for  the  first  time,  a systematic  way  to 
predict  how  oil  spreads  in  calm  water. 

Longuet-Hlgglns  (1969a,  1969b,  1973,  1974)  published  a series  of 
papers  on  the  mect.-’uisms  of  wave  breaking.  He  proposed  g.  nonlinear 
mechanism  for  the  generation  of  sea  waves  and  a model  for  flow  separation 
at  the  free  surface. 

Plate,  et  al.  (1969)  and  Plate  (1970)  measured  the  drift  velocity 
at  the  water  surface  for  two  different  wind  friction  velocities.  They 
observed  that,  over  short  fetches,  the  drift  velocity  increases  according 
to  a power  law,  while,  at  fetches  exceeding  one  meter,  the  drift  velocity 
maintains  a constant  ratio  to  the  wind  velocity. 

Chang,  et  al.  (1970,  1972)  developed  a self-adjusting  probe  posi- 
tioner (wave  follower)  to  position  a crossed-film  sensor  above  the 
instantaneous  water  surface  and,  thus,  were  able  to  measure  the  wind 
velocity  field  as  near  the  instantaneous  water  surface  as  possible. 

Jobson  and  Sayre  (1970a,  1970b)  conducted  a series  of  experiments 
in  a large  open  channel  by  using  a continuous  line  of  dissolved  dye 
across  the  width  of  the  channel.  They  measurod  the  vertical  diffusion 
coefficient. 
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Murray,  et  al.  (1970)  and  Murray  (1972,  1975)  discussed  oceanographic 
observations  near  the  Chevron  spilling  well  off  the  Mississippi  Delta 
and  the  relative  roles  of  the  various  physical  factors  in  the  behavior 
of  oil  slicks.  They  evaluated  the  role  of  large-scale  turbulence  and 
the  effects  of  surface  tension  in  the  control  of  oil  diffusion.  They 
showed  that  the  size  and  shape  of  a given  oil  slick  can  be  predicted 
with  knowledge  of  the  surface  current,  eddy  diffuslvity,  and  the  rate  of 
oil  discharge.  Murray  determined  the  relative  effect  of  local  winds  and 
nesr-aurface  currents  on  determining  the  movement  of  oil  slicks  in 
coastal  and  shelf  waters. 

Prych  (1970)  examined  the  effect  of  density  differences  on  lateral 
mixing  in  open  channels. 

Schwarczberg  (1970,  1971)  Investigated  the  movement  and  spreading 
of  oil  spills  on  water.  Spreading  rates  for  small  spills  were  measured 
and  correlated  with  spill  volume,  oil  density,  and  water  viscosity. 

Field  data  Indicated  that  these  aspects  of  small  spills  cannot  be  scaled 
up  to  large  spill  volumes. 

Spencer  and  Soranson  (1970)  analytically  and  experimentally  studied 
the  setup  cf  oils  having  various  densities  and  viscosities  on  a water 
surface  when  subjected  to  wind.  They  showed  that  the  setup  of  oil  is  a 
function  of  wind  stress  and  wind-generated  wave  action. 

Stewart  (1970)  measured  the  velocity  field  over  deep-water  waves. 

His  measured  spectra  showed  the  presence  of  peaks  caused  by  water  waves. 

Banke  and  Smith  (1971)  presented  spectra  of  all  three  wind  components 
and  co-spectra  for  the  Reynolds  stress  and  sensible  hast  flux  over  water 
in  the  Beaufort  Sea,  in  which  wave  effects  wars  present. 

McNeely  (1971)  described  an  apparatus  for  measuring  oil  slick  dis- 
persion. 

Beebe  (1972)  measured  the  turbulent  flow  over  a wavs  boundary. 

Boardman  and  Dawson  (1972)  presented  a new  process  for  dispsreing 
an  oil  slick. 

Fannelop  and  Waldman  (1972)  considered  the  spread  of  an  oil  slick 
over  calm  water  from  a theoretical  viewpoint.  They  presented  a general 
picture  of  the  slick  development. 
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Hoult  (1972)  wrote,  from  an  analytical  viewpoint,  an  excellent 
review  article  on  oil  spreading  on  the  sea. 

Kondo,  et  al.  (1972)  presented  a spectrum  for  longitudinal  fluctua- 
tions near  the  air-sea  interface  which  indicated  that  an  influence  of 
the  sea  swell  was  present. 

Manton  (1972)  modified  Miles's  wave  theory  to  include  the  effects 
of  the  turbulence  present  in  the  airflow.  A particularly  interesting 
result  was  that  the  pressure  leads  the  water  surface  elevation  in  space, 
and  the  phase  difference  remains  almost  constant  at  75  degrees- 

Rldgway  (1972^  presented  a mathematical  analysis  for  predicting  the 
direction  the  surface  oil  drifts  with  the  wind  and  the  tide. 

Shemdin  (1972)  measured  the  wind-induced  current  and  the  phase 
speed  of  wind  waves. 

Su  (1972)  described  a theoretical  and  experimental  investigation  of 
the  vertical  and  lateral  mixing  of  tracer  fluids  in  a straight,  wide 
open  channel,  with  and  without  wind  on  the  water  surface.  Vertical 
mixing  rates  significantly  increased  when  waves  appeared  on  the  water 
surface. 

Waldman,  et  al.  (1972,  1973)  presented  an  analysis  of  the  spreading 
and  transport  of  oil  slicks  on  the  open  ocean.  They  assessed  the  role 
of  turbulent  diffusion  ss  a transport  mechanism  and  presented  conclusions 
regarding  the  applicability  of  the  various  analytical  models. 

Pottinger  and  Reisbig  (1973)  conducted  laboratory  experiments  to 
determine  the  effect  of  waves  on  oil  spill  drift.  The  data  showed  that 
the  drift  velocity  was  greater  than  that  predicted  by  Stoke's  mass 
transport  theory. 

Reisbig  (1973)  studied  the  effects  of  the  wind  and  waves  on  oil 
spill  drift  so  that  he  could  evaluate  the  importance  of  coupled  wind 
and  wave  drift  mechanisms  in  computer  simulations  of  the  total  drift. 

Reisbig.  et  al.  (1973)  conducted  laboratory  experiments  to  determine 
the  surface  drift  speeds  of  an  oil  lens  which  was  subjected  to  the 
separated  and  coupled  effects  of  wind  waves.  The  results  show  that  the 
wind  drift  and  wave  drift  mechanisms  are  not  simply  additive  over  all 
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regimes  of  wind  speed.  Ac  low  wind  speeds  Che  wave  drift  provides  an 
augmentation  to  the  wind  drift.  At  higher  wind  speeds,  however,  the 
waves  cause  a net  decrease  In  the  coupled  drift  velocity. 

Sturm  and  Sorrell  (1973)  developed  an  optical  system  to  measure 
wave  slope.  The  measured  slope  was  integrated  to  obtain  wave  height. 

Wave  height  results  obtained  from  this  optical  system  were  compared  with 
those  from  a capacitance  wave  height.  The  meniscus  effect  on  the  perfor- 
mance of  the  capacitance  wave  probe  was  determined. 

Veenhuizen  and  Karaki  (1973)  conducted  an  experimental  study  of 
turbulent  airflow  over  water  waves.  They  found  an  organized  wave  motion 
corresponding  to  the  water  wave  in  the  turbulent  airflow.  This  wave- 
induced  motion  was  apparent  up  to  four  or  five  wave  heights  above  the 
mean  water  level. 

Young,  et  al.  (1973)  measured  the  turbulent  heat,  mass  and  momentum 
transfer  mechanisms  in  the  air  boundary-layer  over  wind  waves. 

Banner  and  Phillips  (1974)  investigated  the  incipient  breaking  of 
small-scale  waves.  They  presented  a heuristic  picture  of  the  t>  ’ent 
motion  in  a breaking  wave. 

Kinner , et  al.  (1974)  studied  spoil  disposal  and  oil  dispersion 
experimentally  with  dye. 

Kullenberg  and  Talbot  (1974)  edited  31  papers  on  near-shore  mixing 
caused  by  waves  and  currents,  dispersion  over  the  continental  shelf, 
wind  effects,  oil  slicks,  and  many  other  topics  related  to  the  physical 
processes  responsible  for  the  dispersal  of  pollutants  in  the  sea. 

Lissauer  (1974,  1975)  projected  the  movement  of  oil  slicks  and 
their  impact  iocations  along  the  shorelines  of  different  states  from 
three  potential  deep-water  port  sites.  He  used  average  monthly  wind 
speeds  and  directions,  and  average  monthly  current  patterns  to  predict 
oil  slick  movement. 

Thorton  and  Krapohl  (1974)  measured  wave  height  and  water  particle 
velocities  in  ocean  waves.  The  measured  wave-induced  velocities  were 
comparable  to  those  calculated  by  linear  wave  theory.  At  low  frequencies, 
the  motion  was  almost  totally  wave  induced.  At  higher  frequencies, 
however,  it  was  apparent  that  the  motion  was  primarily  turbulence. 
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WanR  and  Hwang  (1974)  developed  a computer  model  for  simulating  oil 
spreading  and  transport. 

Klemas,  et  al.  (1975a,  1975b,  1976)  conducted  oil-slick  tracking 
experiments  to  verify  a predictive  model  for  oil  dispersion  and  movement. 
They  found  that,  during  certain  parts  of  the  total  cycle,  oil  slicks 
tend  to  line  up  along  estuarine  and  coastal  fronts  or  boundaries. 

Leibovich  (1975a,  1975b)  developed  a statistical  model  to  Investi- 
gate the  effect  of  turbulence  on  removing  oil  from  the  sea  surface.  He 
emphasized  the  importance  of  vertical  turbulent  transport  to  oil  spill 
containment  and  cleanup  operations.  He  also  considered,  theoretically, 
the  performance  of  pil  spill  control  equipment  in  calm  and  rough  water. 
His  analysis  explains  the  observed  thickening  of  oil  in  the  crests  of 
wave 8 . 

Lin  and  Veenhul2en  (1975)  studied  the  effect  of  fetch  on  wind  waves 
by  correlating  experimental  data  obtained  in  the  laboratory  and  in  the 
field.  They  obtained  empirical  formulas  relating  the  air  boundary-layer 
thickness,  the  wind  stress  exerted  on  the  water  surface,  and  the  wave 
height  to  the  wind  speed  and  the  fetch. 

Phillips  and  Banner  (1975)  theoretically  estimated  the  reduction  in 
the  maximum  amplitude  that  short  waves,  superimposed  on  longer  ones, 
attain  when  they  are  at  the  point  of  incipient  breaking. 

Van  Dorn  and  Pazan  (1975)  reported  using  hot-film  probes  success- 
fully near  a breaking  wave. 

Veenhuizen,  et  al.  (1975)  developed  a wind-wave  towing-tank  system 
equipped  with  a servo-controlled  wave  maker  for  generating  random  sur- 
face waves,  a tripping  mechanism  for  generating  a thick  air  turbulent 
boundary-layer  within  a short  distance,  and  an  adjustable  ceiling  for 
controlling  the  pressure  gradient  along  the  tar  k.  This  system  has  been 
used  to  generate  quasi-homogeneous  wind  waves  that  have  a reduced  fetch 
effect. 

Yu  and  Lin  (1975)  developed  a technique  for  generating  a thick, 
turbulent,  air  boundary-layer,  within  a short  fetch,  in  a wind-wave 
tank. 


68 


Danner  and  Melville  (1976)  conducted  laboratory  experiments  to 
investigate  the  airflow  separation  over  a mobile  air-water  interface. 

It  was  argued  that,  in  a frame  of  reference  in  which  the  interfacial 
boundary  assumes  a steady  shape,  the  occurrence  of  separation  requires  a 
etagnation  point  on  the  interface.  They  made  flow  visualization  and 
pressure  measurements  to  demonstrate  air  flow  separation. 

Robben,  et  al.  (1976)  developed  a remote  sensing  technique  for 
measuring  concentrations  of  fluids.  They  us  ’ an  Argon-ion-laser  (1- 
watt)  to  excite  a fluorescent  dye,  Rhodamim.  ' and  measured  the 
fluorescence  intensity  as  an  indication  of  the  dye  concentration.  Dye 
concentrations  as  low  as  10 ^ g/llter  can  be  measured  with  this  technique. 

Stewart  (1976)  studied,  for  application  .o  the  oil  spill  problem, 
the  combined  effects  of  interfacial  elasticity  and  multiple  fluid  layer 
geometties  upon  the  first-order  properties  of  surface  waves. 

Takeuehi  (.1976)  was  the  first  fo  use  a laser  Doppler  veloclmeter 
to  mea«"r‘e  the  turbulence  in  the  water  region  below  the  troughs  of  wind 
waves . 

Thorton,  et  al.  (1976)  measured  waves,  vertical  and  horizontal  water 
particle  velocities  in  spilling,  plunging  and  surfing  breakers  at  sandy 
beaches  in  the  vicinity  cf  Monterey,  Calif.  The  breaking  waves  are 
highly  nonlinear.  Strong  harmonics  appear  in  the  9pectra  and  cross 
spectra  of  wave  and  velocities.  The  strength  of  harmonics  is  related  to 
the  beach  steepness,  wave  height  and  period. 

Donelan  (1977)  measured  the  momentum  transfer  under  breaking  waves 
in  both  laboratory  and  field.  The  momentum  transfer  process  was  quite 
intermittent.  The  immediate  mixing  of  surface  waters  due  to  turbulence 
produced  by  a white  cap  penetrates  to  a depth  of  the  order  of  the  wave 
height . 

Liu,  et  al.  (1977)  applied  the  laser  dye-fluorescence  technique 
to  measure  concentrations  in  a turbulent  Jet. 
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Long  and  Hwang  (1977)  measured  wind  wave  slope  spectra  in  the 
capillary-gravity  range.  Their  results  indicate  that  an  equilibrium 
level  iv  exist  for  each  wind  slope  and  that  it  increases  in  intensity 
with  i sing  wind. 

&lL  (1977)  performed  a theoretical  study  to  determine  the  sea  state 

limit  for  the  survival  of  oil  slicks  on  the  ocean. 

Thornton  (1977)  rederived  the  saturation  range  in  the  frequency 

spectrum  of  wind-generated  waves.  Using  similarity  arguments,  a general 

2 -5 

formulation  was  rederived  giving  Phillips'  wave  energy  spectrum  ~g  f 
in  deep  water.  In  shallow  water,  he  obtained  wave  energy  spectrum 
~ghf  ^ suggested  by  Kitalgordskli  et  al.  (1975).  The  velocity  Bpectrum 
in  the  saturation  range  is  derived  using  linear  wave  theory  as  a spectral 
transfer  function.  In  both  deep  and  shallow  water  he  obtained  the  hori- 
zontal velocity  spectrum  -f  ^ . 
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APPENDIX  B.  COMPARISON  OF  MEASUREMENTS  BY  LDV  AND  HOT-FILM  PROBES 


In  Sections  4.5  and  4.6,  ve  have  presented  the  measurements  of  mean 
velocity  and  rms  velocity  fluctuations  in  wind  waves  using  an  array  of 
hot-film  probes.  Before  we  proceeded  to  use  the  LDV  to  measure  the 
velocity  field  in  mechanically-generated  waves  with  wind,  we  tested  the 
LDV  performance  in  wind  waves,  and  made  a comparison  of  measurements  by 
the  LDV  and  hot-film  probes.  The  experimental  conditions  for  the  com- 
parison measurement  were  the  same  as  described  in  Section  4 except  that 
the  LDV  measurement  is  positioned  at  x " 4.75  m and  the  measurement 
points  were  20  cm  away  from  the  side  wall  while  the  hot-film  probes  were 
positioned  at  x ■ 5 m along  the  centerline  of  the  tank  which  is  58  cm 
from  the  side  wall.  The  LDV  was  positioned  at  x * 4.75  m because  at 
x ■ 5 m there  was  a vertical  column  blocking  the  laser  light  from  the 
optics  of  the  LDV.  We  set  the  measurement  point  for  the  LDV  20  cm  away 
from  the  side  wall  because  there  is  a limitation  in  the  scattering 
signal  required  to  obtain  good  record  of  velocity  by  the  LDV.  In  other 
words,  if  the  focal  distance  is  too  large,  the  scattering  becomes  so 
weak  that  the  LDV  measurement  is  difficult.  Further  optimization  of  the 
present  system  is  feasible  to  increase  the  focal  distance. 

Figures  72  and  73  present  the  comparisons  for  mean  velocity  and 
v is  longitudinal  velocity  fluctuations.  The  agreement  between  these  two 
measurements  is  very  good  in  view  of  the  fact  that  they  were  indepen- 
dently made  with  two  different  instruments  at  different  times.  This 
comparison  establishes  our  confidence  in  both  instruments,  the  LDV  and 
the  hot-film  probe. 
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APPENDIX  C.  EFFECTS  OF  AN  OIL  SLICK  ON  WIND  WAVES 

Experiments  were  conducted  to  investigate  the  effects  of  an  oil 

slick  on  wind  waves  for  various  wind  speeds  and  at  several  fetches. 

Diesel  oil  No.  2 was  used  In  the  experiments.  It  was  fed  onto  the  water 

surface  from  the  upstream  end  of  the  wave  tank,  and  the  discharge  was 

0.35  llter/aec.  The  specific  gravity  of  the  oil  was  0.84  and  the 

2 

kinematic  viscosity  was  0.19  cm  /sec. 

The  photodiode  wave  gauge  was  used  to  measure  the  vertical  displace- 
ment of  the  water  surface  with  and  without  an  oil  slick  present.  The 
gauge  was  positioned  along  the  center  line  of  the  wave  tank.  In  case  an 
oil  slick  was  present,  the  intensity  of  the  oil  flourescence  was  strong 
enough  that  the  interface  between  the  air  and  the  oil  could  be  detected 
with  the  photodiode  wave  gauge.  In  all  the  experiments  reported  here, 
fluorescein  disodium  dye  was  added  to  the  water  so  that  the  interface  of 
the  water  and  the  air  could  be  detected  when  the  oil  slick  was  absent. 

For  each  experimental  run,  we  first  took  data  for  tne  case  without  an 
oil  slick.  Then  an  oil  slick  was  fed  onto  the  water  surface  to  investi- 
gate Its  effects  on  wind  waves. 

Figure  74  shows  a record  of  the  vertical  displacement  of  the  water 
or  oil  surface  at  a fetch  x • 5 m for  wind  speed  - 10  m/sec.  Before 
the  oil  slick  reaches  the  measurement  point,  the  record  shows  the  time 
variation  of  the  vertical  displacment  of  the  water  surface.  As  the  oil 
alick  reaches  the  measurement  point,  the  mean  level  of  the  air-oil 
Interface  Increases  and  the  fluctuating  vertical  displacement  decreases. 
The  transient  stage  lasts  for  about  three  to  four  minutes  and  a steady 
stage  reestablishes  four  minutes  afterwards.  The  increase  of  the  mean 
level  is  partly  due  to  the  thickness  of  the  oil  slick  and  partly  due  tc 
the  fact  that  the  setup  of  the  water  surface  Is  altered  by  the  presence 
of  the  oil  slick.  After  the  steady  stage  is  reached,  another  set  of 
data  is  recorded  for  the  case  when  the  oil  slick  is  present.  This  set 
of  data  could,  therfore,  be  compared  with  data  obtained  without  an  oil 
slick  present  in  order  to  investigate  the  effects  of  oil  slick  on  wind 
waves , 
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C.l  Wave  Profiles 


Figure  75  shows  wind-wave  profiles  without  an  oil  slick  present, 
which  were  measured  at  fetch  x * 4 m for  three  wind  speeds,  U,  - 4-,  7-, 
and  10  m/sec.  The  wind-wave  profiles  are  plotted  as  the  time  variation 
of  the  vertical  displacement  of  the  air-water  interface.  The  amplitude 
and  period  of  dominant  waves  increase  as  wind  speed  increases  and  wave 
packets  are  observed.  These  results  are  consistent  with  those  presented 
in  Section  4.2. 

Figure  76  shows  the  wind  wave  profiles  when  an  oil  slick  is 
present.  Comparison  between  Figures  75  and  76  indicates  that  the 
wind  waves  are  significantly  damped  by  the  oil  slick.  In  Figure  76(c), 
small  scale  fluctuations  appear  which  indicate  the  break  up  of  the  oil 
slick.  Visually  we  have  seen  the  formation  of  oil  lenses  floating  on 
the  water  surface  for  ■ 10  m/sec. 

C. 2 RMS  Wave  Height 

The  rras  values  of  the  wave  displacement  are  plotted  versus  wind 
speed  in  Figure  77.  When  the  oil  slick  is  absent,  the  rms  wave  height 
increases  with  wind  speed  and  fetch,  which  is  consistent  with  the  results 
presented  in  Section  4.2.  For  all  wind  speeds,  the  rms  wave  height  for 
the  case  with  an  oil  slick  present  is  smaller  than  that  for  the  case 
without  an  oil  slick  present.  Thus,  the  damping  effects  of  an  oil  slick 
on  wind  waves  are  evident.  When  an  oil  slick  is  present,  however,  the 
rms  wave  height  increases  when  going  from  4-  to  7 m/sec  wind,  but 
decreases  when  going  from  7-  to  10  m/sec  wind.  It  is  worth  mentioning 
that  for  the  case  of  wind  speed  at  10  m/sec,  oil  lenses  formed  on  the 
water  surface.  We  speculate  that  these  oil  lenses  could  transfer  energy 
from  air  to  water  more  effectively  than  a continuous  sheet  of  oil  slick 
because  the  presence  of  oil  lenses  makes  the  water  surface  rougher,  and 
the  waves  could  be  more  effectively  damped  by  them.  Figure  78  shows  the 
variation  of  rms  wave  height  with  fetch.  For  the  case  without  an  oil 
slick,  the  rms  values  increase  with  fetch,  which  is  consistent  with  the 
results  obtained  in  Section  2.  When  an  oil  slick  is  present,  the  rms 
values  do  not  vary  appreciably  with  fetch. 
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C. 3 Wave  Energy  Spectrum 

Figure  79  presents  the  wave  energy  spectrum  S(f)  measured  with  the 
photodiode  wave  gauge  located  at  fetch  x *>  3.2  m for  wind  speed  • 10 
m/sec  without  an  oil  slick  present.  Two  distinct  frequency  regimes  are 
observed,  namely,  a gravity  wave  regime  and  a capillary  wave  regime. 

The  gravity  waves  have  a peak  spectrum  at  a frequency  of  about  4.7  Hz. 

In  the  capillary  wave  regime,  the  wave  energy  spectrum  follows 

S(f)  - $ y2/3  f“7/3  , 7.5  Hz  < f < 100  Hz  , 

where  y is  the  surface  tension  per  unit  density  and  has  a value  of 

3 2 

70  cm  /sec  for  water,  and  S is  a numerical  constant  having  a value  of 

_ 3 

about  7.2  x 10  for  this  particular  measurement.  This  -7/3  power  law 
was  obtained  based  on  dimensional  arguments  (Phillips,  1966),  and  it  was 
suggested  originally  by  Hicks  (See  Phillips,  1958).  To  the  knowledge  of 
the  authors,  the  wave  energy  spectrum  presented  in  Figure  79  is  the 
first  experimental  verification  of  this  power  law.  The  measurement  of 
the  wave  energy  spectrum  at  high  frequencies  is  very  important  in  under- 
standing the  characteristics  of  capillary  waves.  Since  these  capillary 
waves  act  as  roughness  elements  to  wind,  their  characteristics  would 
Influence  the  wind  stress  exerted  on  the  water  surface  and,  therefore,  the 
surface  drift  current.  In  Figure  79,  there  are  two  spectra  calculated 
from  two  Independent  records.  The  shaded  area  provides  an  estimate  of 
statistical  error  In  measuring  the  wave  energy  spectrum. 

The  wave  energy  spectra  measured  at  a fetch,  x » 4u  for  three  wind 
speeds,  ■ 4-,  7-,  and  10  m/sec  without  an  oil  slick  present  are  plot- 
ted in  Figure  80.  Each  spectrum  has  two  distinct  wave  regimes  as 
aforementioned.  The  value  of  8 increases  with  wind  speed.  The  dominant 
frequency  decreases  with  wind  speed.  Also,  the  frequency  at  which  two 
distinct  wave  regimes  Join  decreases  with  wind  speed.  Thus  the  measure- 
ments show  that  as  wind  speed  increases  the  wave  energy  in  both  gravity 
and  capillary  wave  regimes  increases.  Long  and  Huang  (1977)  measured 
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wind-wave  slope  spectra  in  the  capillary  and  gravity  ranges.  Their 
results  indicate  that  an  equilibrium  level  may  exist  for  each  wind 
speed,  and  that  it  increases  in  intensity  with  increasing  wind  speed. 

The  present  results  are  therefore  consistent  with  theirs.  It  would  be 
interesting  to  investigate  whether  an  asymptotic  nature  of  the  wave 
energy  spectrum  exists  in  the  capillary  wave  regime. 

The  effects  of  an  oil  slick  on  wind  waves  were  also  investigated 
from  measurements  of  wave  energy  spectra.  Figure  81  shows  the  wave 
energy  spectra  measured  at  fetch  x * 4 m for  wind  speed  ■ 4 m/stec, 
with  and  without  an  oil  slick  present.  The  oil  slick  influences  the 
wave  spectrum  over  all  wave  frequencies.  In  particular,  the  dominant 
wave  frequency  for  the  case  of  an  oil  slick  being  present  is  higher  than 
that  for  the  case  without  an  oil  slick.  This  result  indicates  that  the 
onset  of  the  wind  waves  is  delayed  by  the  presence  of  an  oil  slick. 

Figure  82  shows  the  wave  energy  spectra  measured  at  x • 4 m for 
Uro  a 7 m/sec.  When  an  oil  slick  is  present,  the  energy  of  the  dominant 
waves  decreases,  but  the  dominant  wave  frequency  remains  unchanged.  At 
higher  frequencies,  the  effect  of  an  oil  slick  is  not  as  appreciable  for 
- 4 m/sec.  Comparison  between  Figures  81  and  82  shows  that  the 
effect  of  an  oil  slick  depends  on  wind  speed.  Figure  83  shows  a com- 
parison of  wave  energy  spectra  measured  with  and  without  an  oil  slick 
present  for  wind  speed  U ■ 10  m/sec.  The  spectral  density  at  the 
dominant  wave  frequency  decreases  significantly  when  the  oil  slick  is 
present.  In  the  frequency  range  6 Hz  < f < 20  Hz  , the  spectral  density 
appears  flat  and  it  is  possibly  contributed  by  the  presence  of  oil 
lenses,  lor  frequencies  f > 20  Hz  , the  oil  slick  effect  is  not 
appreciable.  From  Figures  82  and  83,  we  conclude  the  following  for  the 
present  experimental  conditions.  For  wind  speed  > 7 m/sec,  the 
effect  of  an  oil  slick  on  the  generation  of  capillary  waves  is  not 
significant.  Tu?ae  capillary  waves  most  likely  are  generated  locally  by 
wind  and  governed  by  surface  tension.  Therefore,  the  evolution  of 
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capillary  waves  is  expected  to  be  different  from  that  of  gravity  waves. 
More  research  is  needed  to  explore  this  evolution. 
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APPENDIX  D.  EFFECTS  OF  AN  OIL  SLICK  ON  MECHANICALLY-GENERATED  WAVES  WITH 

WIND 

Experiments  were  conducted  to  Investigate  the  effects  of  an  oil 
slick  on  mechanically-generated  waves  with  wind.  The  wave  maker  was 
oscillated  at  a frequency  of  1.6  Hz,  and  three  different  strokes, 

Kq  ■ 2.54-,  5.08-,  and  7.62  cm,  were  used.  The  Wind  speed  was  10  m/sec 
except  for  a few  cases  in  which  wind  speeds  of  4-  and  7 m/sec  were  used. 
Measurements  of  wave  displacement  were  made  with  the  photodiode  wave 
gauge  at  four  fetches,  x ■ 2-,  3.2-,  4-  and  5 m.  For  both  the  cases 
with  and  without  an  oil  slick,  the  experimental  procedure  was  the  same 
as  described  in  Appendix  C.  Experimental  results  of  wave  profiles,  rms 
wave  heights  and  wave  energy  spectra  were  obtained  to  examine  the 
effects  of  oil  slick. 

D. 1 Wave  Profiles 

Figure  84  is  the  wave  profile  measured  at  x ■ 5 m for  Hq  ■»  2.54-, 
5.08-,  and  7.62  cm  and  ■ 10  m/sec.  In  this  figure,  the  apparent 
propagation  direction  of  the  waves  is  from  the  right  to  the  left.  The 
vertical  scale  is  exaggerated  to  show  the  details  of  the  profiles.  All 
three  profiles  have  the  dominant  wave  frequency  corresponding  to  the 
oscillating  frequency  of  the  wave  maker,  i.e.,  f ■ 1.6  Hz.  Also, 
superimposed  on  the  dominant  waves  are  higher  harmonic  wave  components 
and  high  frequency  wave  components  generated  by  wind.  When  the  stroke 
of  the  wave  maker  Increases,  the  wave  steepness  increases,  and,  as  a 
result,  the  waves  become  unstable  and  begin  to  break.  For  Hq  - 5.08  cm, 
wave  breaking  appeared  as  indicated  by  the  small-scale  waves  resting  on 
wave  crests,  and  for  Hq  ■ 7:62  cm  the  breaking  is  more  violent. 

In  Figures  85  and  86,  we  replot  the  wave  profile  to  scale  approxi- 
mately. Figure  85  shows  the  case  for  ■ 2.54  cm,  and  the  wave  profile 
does  not  appreciably  deviate  from  a sinusoidal  profile.  Small-scale  or 
high  frequency  wave  components  are  generated  mostly  by  wind  and  their 
profile  is  shown  in  the  enlarged  plots.  Wave  breaking  does  not  occur 
yet.  Figure  86  shows  the  wave  profile  for  Hq  ■ 7.62  cm.  The  wave 
profiles  are  highly  asymmetric  because  the  wave  slope  in  the  downstream 
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(downwind)  side  is  much  larger  chan  that  in  the  upstream  side,  and  the 
wave  slope  at  the  wave  crest  is  steeper  than  that  at  the  wave  trough. 

Two  enlarged  wave  profiles  at  wave  crests  are  re-plotted,  and  they 
represent  two  realizations  of  wave  breaking.  These  plots  clearly  show 
that  the  breaking  process  is  quite  random  as  these  two  realizations  do 
not  resemble  each  other  at  all. 

Figure  87  shows  the  wave  profiles  measured  at  x = 5m  for  Ho  = 2.54-, 
5.08-,  and  7.62  cm  and  * 10  m/sec  when  an  oil  slick  is  present. 

Direct  comparions  between  Figures  87  and  84  reveal  the  effects  of  an  oil 
slick  on  waves.  For  Hq  • 2.54  cm,  the  overall  wave  amplitude  decreases, 
and  higher  harmonic  wave  components  disappear  when  an  oil  slick  is 
present.  For  Hq  “ 5.08  cm,  the  presence  of  an  oil  slick  makes  the  wave 
slopes  at  crests  steeper.  The  wave  breaking  is  somehow  inhibited  by 
the  oil  slick.  As  shown  in  Figure  87  (b) , small-scale  wave  components 
hardly  exist  when  comparing  with  Figure  84  (b).  For  Hq  ■ 7.62  cm,  wave 
breaking  occurs  at  wave  crests.  Comparison  between  Figures  84c  and  87c 
shows  that  the  wave  breaking  is  so  violent  that  it  is  not  affected 
significantly  by  the  oil  slick. 

D.2  RMS  Wave  Height 

The  rtns  wave  heights  measured  with  and  without  an  oil  slick  are 
shown  in  Figure  88.  The  data  plotted  in  the  figure  include  three  wave 
conditions,  H « 2.54-,  5.08-  and  7.62  cm,  measured  at  four  fetches, 
x = 2-,  3.2-,  4-  and  5 m.  Each  rms  wave  height  was  calculated  from  a 
record  of  about  60  seconds. 

With  the  oil  slick  present,  the  wave  energy  still  grows  with  fetch, 
but  at  a slower  rate  than  that  without  the  oil  slick  present  for  r 
2.54  cm.  For  Hq  *•  5.08  cm,  the  effect  of  the  oil  slick  on  the  wave 
growth  is  not  significant.  For  «*  7,62  cm,  however,  the  wave  energy 
decreases  as  the  fetch  Increases  from  x « 2 m to  x m 3. 2 m,  and  then 
increases  slowly  for  the  fetch  varying  from  x “ 3.2  m to  x * 5 m.  Flow 
visualization  results  show  that  at  small  fetches,  say  x ■ 2 m,  the  oil 
slick  breaks  into  oil  lenses  by  the  wave  motion  and  turbulence.  The  oil 
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lenses  which  float  on  the  water  surface  could  damp  waves  very  effectively 
(see  Appendix  C.2).  Further  downstream,  the  oil  lenses  are  broken  into 
oil  droplets  which  in  turn  disperse  into  the  water  column  by  the  wave 
motion  and  turbulence.  Thus,  the  oil  concentration  near  the  water 
surface  becomes  diluted,  and  the  effect  of  the  oil  slick  is  no  longer  as 
strong  as  at  small  fetches.  The  measurements  of  the  rms  wave  height  are 
plotted  versus  wind  speed  at  x * 5 m for  Hq  - 7.62  cm  in  Figure  89. 

At  U ■ 4 and  7 m/sec  the  wave  energy  measured  with  the  oil  slick  is 
smaller  than  that  measured  without  the  oil  slick.  But,  at  IJ^  ■ 10  m/sec, 
the  effect  of  oil  slick  is  not  measureable,  because  under  this  condition 
most  of  the  oil  slick  is  transported  vertically  below  the  water  surface. 

D. 3 Wave  Energy  Spectrum 

In  this  section,  the  wave  energy  spectra  measured  with  and  without 
the  oil  slick  are  compared.  Attempts  are  made  to  clarify  further  the 
effects  of  an  oil  slick  on  mechanically-generated  waves  with  wind. 

In  Figure  90,  comparison  of  two  spectra  measured  with  and  without 

the  oil  slick  for  H * 2. 54  cm  at  x ® 5 m is  shown.  The  oil  slick  damps 

o 

cut  considerably  the  wave  components  with  f < 20  Hz  . 

For  = 5.08  cm  (Figure  91),  the  spectral  density  is  higher  at 

f ■ f , 2f  and  3f  with  the  oil  slick  than  without.  This  confirms  the 
o o o 

steepening  of  the  wave  crests,  which  is  observed  in  Figure  87,  due  to 

the  reduction  in  surface  tension  when  the  oil  slick  is  present.  As  a 

result,  the  rms  wave  height  is  enhanced  slightly  rather  than  damped. 

With  further  Increase  of  H 7.62  cm,  no  noticeable  difference 

o 

between  the  two  spectra  is  observed  in  Figure  92.  Note  that  intensive 
breaking  is  In  progress  in  both  cases.  The  oil  slick  is  broken  into 
small  droplets  and  a great  number  of  the  droplets  are  transported  into 
the  water  column.  As  a result,  the  effect  of  the  oil  slick  in  wave 
damping  is  insignificant.  This  result  thews  that  for  « 7.62  cm,  and 
■ 10  m/sec,  the  oil  slick  could  be  treated  as  a passive  medium  as  it 
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barely  affects  the  wave  structure.  In  other  words,  the  effects  of  oil 
viscosity  and  surface  tension  are  unlikely  to  be  important  in  the 
present  study. 
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TABLE  1.  -SUMMARY  OF  OIL  PROPERTIES,  ARCO  DIESEL 
(ASTM  Grade  2D) 


VISCOSITY 

TEMP. 

SPECIFIC 

DATE 

(CPS) 

<°C) 

GRAVITY 

12-12-77 

18.8 

20.2 

0.838 

12-14-77 

18.9 

20.1 

0.840 

12-15-77 

- 

19.8 

0.840 

12-16-77 

- 

- 

0.840 

12-19-77 

19.6 

19.3 

0.842 

12-21-77 

19.0 

- 

0.842 

12-21-77 

19.5 

18.9 

0.844 

12-22-77 

19.6 

19.4 

0.844 

1-4-78 

20.2 

18.4 

0.845 

1-5-78 

20.0 

19.4 

0.845 

1-6-78 

20.2 

19.5 

0.846 

1-9-78 

20.4 

19.7 

0.845 

1-10-78 

20.5 

17.9 

0.846 

1-11-78 

20.3 

19.4 

0.846 

1-13-78 

20.9 

17.5 

0.848 

1-16-78 

20.2 

19.7 

0.846 

1-17-78 

20.4 

19.3 

0.846 

1-18-78 

20.3 

18.4 

0.848 

1-19-78 

18.5 

19.0 

0.840 

1-19-78 

- 

- 

0.847 

1-20-78 

20.0 

19.7 

0.846 

1-23-78 

20.4 

19.0 

0.844 

1-24-78 

18.8 

19. 1 

0.844 

1-25-78 

19.0 

19.0 

0.842 

1-26-78 

18.9 

19.  1 

0.840 

1-27-78 

18.8 

17.5 

0.840 

1-30-78 

18.9 

18.3 

0.841 

1-31-78 

18.5 

18.9 

0-840 

2-1-78 

18.8 

18.2 

0.842 

2-2-78 

19.0 

19. 1 

0.841 

2-3-78 

19.2 

19.  1 

0.841 
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TABLE  2 . -WATER 

SURFACE 

DRIFT  VELOCITY  FOR 

WIND-WAVES 

RUN 

NUMBER 

(tn/sec) 

X 

(m) 

NO.  OP  FLOATS 
OBSERVED 

(cm/«ec) 

STANDARD 

DEVIATION 

(cm/eec) 

OIL-34 

4 

1.75 

11 

12.8 

0.8 

35 

4 

2.25 

6 

11.5 

1.2 

36 

4 

3.21 

10 

13.2 

1.4 

37 

4 

4.00 

8 

13.2 

1.5 

38 

4 

4.00 

8 

13.6 

0.5 

44 

4 

1.75 

11 

11 .9 

1.0 

45 

4 

2.25 

4 

10.2 

0.3 

47 

4 

1.75 

31 

11.6 

1.3 

48 

4 

2.25 

17 

10.3 

1.3 

49 

4 

3.21 

14 

12.2 

1.0 

50 

4 

4.00 

16 

12.2 

i.: 

51 

4 

4.75 

13 

13.3 

0.8 

52 

4 

5.25 

14 

12.7 

1.0 

53 

4 

6.00 

14 

12.9 

1.4 

54 

7 

1.75 

25 

20.8 

1.2 

55 

7 

2.25 

23 

21.8 

1.1 

56 

7 

3.21 

15 

22.9 

1.3 

57 

7 

4.00 

15 

25.0 

1.9 

58 

7 

4.75 

21 

24.5 

2.3 

59 

7 

5.25 

2) 

23.5 

2.1 

60 

7 

6.00 

28 

23.8 

5.8 

61 

10 

1.75 

20 

30.8 

3.1 

62 

10 

2.25 

25 

33.3 

2.8 

63 

10 

3.21 

36 

32.5 

3.5 

64 

10 

4.00 

31 

35.0 

5.3 

65 

10 

4.75 

24 

35.2 

6.8 

66 

10 

5.25 

27 

32.6 

3.1 

67 

10 

6.00 

21 

34.1 

3.4 
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TABLE  3.  -WATER  BOUNDARY-LAYER  THICKNESS  FOR  WIND-WAVES 


RUN 

NUMBER 

U 

oo 

(m/ sec) 

X 

(m) 

u; 

NO.  OF 

OBSERVATIONS 

6 

w 

(cm) 

STANDARD 

DEVIATION 

(cm) 

&§* 

U* 

CD 

OIL-68 

4 

1.7 

1.04 

4 

1.24 

0. 17 

7.6  * 

: 10~3 

71 

4 

2.3 

1.41 

10 

0.75 

0.17 

4.6  * 

: 10~3 
-3 

74 

4 

3.2 

1.97 

10 

0.80 

0.13 

4.91 

x 10 
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a.  WIND  TUNNEL 


b.  WAVE  TANK 


FIGURE  1 . WIND-WAVE  TANK 
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FIGURE  6.  PLIOLITE  PARTICLES  PHOTOGRAPHED  WITH  ELECTRONIC 
MICROSCOPE. 
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FIGURE  8.  COMPARISON  OF  LDV  OUTPUT  WITH  PITOT  TUBE  OUTPUT. 
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FIGURE  9 CALIB 
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FIGURE  14.  WAVE  PEOPLE  MEASURED  WITH  PHOTODIODE  WAVE  GAUGE. 


FIGURE  15.  WAVE  PROFILE  MEASURED 


Ha 


U (cm/sec) 

FIGURE  1 8.  MEAN  VELOCITY  PROFILE  IN  AIR  BOUNDARY-LAYER. 


FIGURE 1 9.  INSTANTANEOUS  WATER  SURFACE  DISPLACEMENT  OF  WIND-WAVES. 


FIGURE  21 . VARIATION  OF  RMS  WAVE  AMPLITUDE  WITH  WIND  SPEED. 


FIGURE  23.  FETCH  EFFECTS  ON  DOMINANT  WAVE  FREQUENCY. 
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FIGURE  25.  FETCH  EFFECTS  ON  WAVE  PHASE  SPEED. 
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FIGURE  26.  VARIATION  OF  DOMINANT  WAVE  LENGTH  WITH  WIND  SPEED. 


x(rn) 

FIGURE  29.  VARIATION  OF  SURFACE  DRIFT  VELOCITY  WITH  FETCH. 


=|=—  Water 
Surface 


(b)  Uoo  = 4 m/sec , x = 3.2  m 

FIGURE  30.  DYED  WATER  BOUNDARY-LAYER. 
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(c)  Uuo  = 7 m/sec , x = 1.7m 


(d)  Uoo  = 7 m/sec  , x = 3.2  m 


FIGURE  30.  DYED  WATER  BOUNDARY-LA'  ER.  ICONT.) 


1 30 


Water 
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FIGURE  30.  OYED  WATER  BOUNDARY-LAYER.  (CONT.) 
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FIGURE  32.  MEAN  VELOCITY  PROFILE  IN  WATER  BOUNDARY-LAYER. 
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FIGURE  34. 
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FIGURE  36.  RMS  OF  LONGITUDINAL  VELOCITY  FLUCTUATIONS  IN  WATER 
BOUNDARY-LAYER. 
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FIGURE  38.  RMS  OF  VERTICAL  VELOCITY  FLUCTUATIONS  IN  WATER  BOUNDARY-LAYER. 
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FIGURE  46.  WATER  SURFACE  PHOTOGRAPH. 
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FIGURE  47.  AIR  BUBBLES  ENTRAINED  INTO  WATER  UNDER  WAVE  CREST. 
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FIGURE  48.  WAVE  ENERGY  SPECTRA  (CAPACITANCE  PROBE). 
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FIGURE  60.  VERTICAL  PROFILES  OF  RMS  LONGITUDINAL  VELOCITY  IN 
BREAKING  AND  NONBREAKING  WAVES. 
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FIGURE  01 . VERTICAL  PROFILES  OF  RMS  VERTICAL  VELOCITY  IN  BREAKING 
AND  NONBREAKING  WAVES. 
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FIGURE  62.  VERTICAL  PROFILES  OF  RMS  VELOCITY  FLUCTUATIONS  IN 
BREAKING  AND  NONBREAKING  WAVES. 
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FIGURE  63.  VERTICAL  PROFILERS  MEAN  SQUAF  vi  )<J! 
DERIVATIVES  IN  BREAKING  AND  NONBREAKING  WAV’  b 
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HGURE  64.  SPECTRUM  OF  LONGITUDINAL  VELOCITY  FLUCTUAT  IONS  IN 
BREAKING  WAVES. 
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FIGURE  65.  SPECTRUM  OF  VERTICAL  VELOCITY  FLUCTUATIONS  IN 
BREAKING  WAVES. 
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FIGURE  67.  SPECTRUM  OF  VERTICAL  VELOCITY  FLUCTUATIONS  IN 
BREAKING  WAVES. 
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FIGURE  70.  VERTICAL  PROFILE  OF  MEAN  OIL  VOLUME  CONCENTRATION. 
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FIGURE  80.  WAVE  ENERGY  SPECTRUM  VARIATION  WITH  WIND  SPEED. 
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FIGURE  83.  WAVE  ENERGY  SPECTRA  MEASURED  WITH  AND 
WITHOUT  OIL. 
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FIGURE  87.  PROFILE  OF  MECHANICAL  GENERATED  WAVES  WITH  WIND 
AND  OIL  STICK.  VERTICAL  SCALE  = 2 cm/div.(  HORIZONTAL 
SCALE  = o.5  sec/div. 
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FIGURE  88  VARIATION  OF  WAVE  HEIGHTS  WITH  FETCH. 


FIGURE  89.  VARIATION  OF  RMS  WAVE  HEIGHTS  WITH  WIND  SPEED. 
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FIGURE  82.  COMPARISON  OF  WAVE  ENERGY  SPECTRA  WITH  AND 
WITHOUT  OIL  SLICK. 
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